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1 . 0 INTRODUCTION 

Depending  upon  its  nature,  an  upstream  disturbance  may 
induce,  or  suppress,  a subsequent  downstream  flow  separation. 

There  are  countless  practical  examples  of  such  situations  (from 
the  golf  ball  to  the  airfoil) . The  present  study  is  devoted  to 
understanding  such  flows  at  transonic  speeds.  The  motivation  of 
this  study  is  clearly  explained  in  the  first  phase  summary  report 
(Ref.  1).  However,  it  is  still  worthwhile  to  reiterate  some 
crucial  points  so  vital  in  understanding  today's  transonic  flights. 

It  was  Pearcey  (Ref.  2)  who  pointed  out,  after  an  extensive 
experimental  study,  that  the  flow  over  a transonic  airfoil  or 
wing  presents  two  kinds  of  distinguished  flow  separation  patterns 
over  the  surface.  One  kind  is  a separation  bubble  which  forms  at 
the  foot  of  a near-normal  shock  wave  adjacent  to  the  surface.  The 
bubble  will  grow  aft  as  either  the  free-stream  Mach  number  or  the 
angle  of  attack  increases.  He  called  this  type  of  flow  separation 
pattern,  Type  A.  The  other  kind  is  a flow  separation  from  the 
trailing  edge,  with  or  without  the  local  shock  induced  separation 
near  the  midchord.  This  type  of  separation  pattern  is  classified 
as  Type  B.  In  order  to  see  more  clearly  these  two  types  of  flow 
field  developments,  Pearcey' s original  concept  has  been  reillus- 
trated as  shown  in  Figure  1.  Whether  it  belongs  to  Type  A or  Type 
B,  the  "upstream  influence"  to  the  downstream  flow  is  very  well 
explained.  The  problem  becomes  serious  if  it  is  Type  B,  because 
the  downstream  flow  is  associated  with  a flow  separation.  As  is 
explained  in  this  report,  the  upstream  disturbance  definitely 
affects  the  downstream  flow  condition  and  thus  extends  its 
influence  to  the  flow  separation. 
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(a-1)  Bubble  initially 
induced 


(a-2)  Rear  separation 
starting 


conjunction  tion  moving 

with  rear  upstream 

separation 


\ ♦ ✓ 


(c)  Separation  from  shock 
to  trailing  edge 

Figure  1.  Schematic  diagram  of  typical  transonic  flow 
inviscid/viscous  interactions  (adapted  from 
Ref.  2) . 
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Recently  Yoshihara  and  Zonnars  (Ref-  3)  reported  that  the 
three-dimensional  transonic  shock-induced  separation  over  a super- 
critical wing-fuselage  configuration  is  very  much  a Type  B pat- 
tern (See  Figure  2).  Whitcomb  (Ref.  4)  indicated  also  that  for 
flow  over  a supercritical  wing  with  a contoured  fuselage  (designed 
by  area  rule),  the  flow  separation  pattern  definitely  belongs  to 
Type  B.  The  surface  oil  flow  trace  showed  that  the  flow  started 
to  separate  from  the  trailing  edge  as  the  angle  of  attack  was 
increased.  This  and  the  work  of  additional  researchers  indicate 
that  the  flow  over  a modern  super- critical  wing-fuselage  configu- 
ration is  of  the  Type  B and  that  the  upstream  disturbance  affects 
the  downstream  separation.  Therefore,  the  present  study  is  very 
timely  in  shedding  some  light  in  understanding  this  vital  problem. 

Whenever  flow  separation  is  influenced  by  some  upstream  dis- 
turbance, a large  interaction  can  ensue.  Both  the  rate  of 
relaxation  from  the  upstream  disturbance  and  the  upstream  propa- 
gation from  the  separation  will  be  modified  due  to  the  interaction 
between  the  two  flow  phenomena.  The  local  flow  taking  place 
between  the  upstream  disturbance  and  the  downstream  separation 
will  also  be  changed.  The  extent  of  this  change  must  evidently 
depend  upon  the  following  parameters: 

1.  Free  stream  conditions --Mach  number,  Reynolds  Number, 
etc . 

2.  The  strength  of  the  upstream  disturbance  which  can  be 
characterized  in  many  different  ways.  In  the  absence 
of  a downstream  separation,  perhaps  the  most  convenient 
parameters  for  the  description  of  the  strength  of  the 
upstream  disturbance  are:  the  magnitude  of  the  maximum 
recovered  pressure,  the  boundary  layer  integral  thick- 
ness, and  the  intensity  of  the  turbulence  at  certain 
locations  downstream  of  the  reattachment  point. 
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Figure  2.  Shock  and  separation  pattern  at  45°  sweep, 
Mw  = 0.95  and  a = 9°  (Ref.  3). 
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3.  The  strength  of  the  downstream  separation  which  can  be 
characterized  by  such  parameters  as  the  magnitude  of 
the  peak  pressure  in  the  separation  region,  the  boundary 
layer  integral  thicknesses,  and  the  intensity  of  the 
turbulence  at  the  separation  point. 

4.  The  distance  between  the  upstream  disturbance  and  the 
downstream  separation.  The  flow  after  the  disturbance 
needs  a certain  relaxation  distance  to  recover. 

5.  The  manner  of  distortion  in  the  pressure  and  the  velocity 
profile  as  well  as  their  recovery  process. 

The  present  interim  report  does  not  provide  conclusive  under- 
standing of  the  influence  of  all  the  parameters  listed  above; 
however,  some  significant  progress  was  obtained.  For  instance, 
the  relaxation  distance  (as  mentioned  in  point  4 above)  for  the 
pressure  was  found  to  be  much  shorter  than  that  for  the  velocity 
profile  and  the  pressure  and  skin  friction  recovered  over  different 
distances.  Once  the  pressure  and  skin  friction  recovered,  they 
reached  higher  values  than  the  no  disturbance  condition.  It  has 
been  observed  by  Whitcomb  in  wind  tunnel  tests  (Ref.  4)  that  a 
constant  pressure  region  for  a small  distance  aft  of  the  shock 
wave  on  a supercritical  airfoil  actually  helps  in  preventing 
trailing  edge  separation  (Figure  3).  His  observations  can  be 
clearly  explained  from  the  findings  reported  herein. 

The  basic  configuration  used  for  the  experiments  conducted 
in  this  study  is  shown  in  Figure  4.  This  figure  shows  the 
nomenclature  used  in  this  report  and  gives  a schematic  drawing 
of  the  pressure  distribution  and  the  flow  over  a shallow-cavity 
model . 
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Figure  3.  Comparison  of  pressure  distribution  on 
conventional  and  supercritical  airfoils. 
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2.0  FACILITIES  AND  EXPERIMENTAL 
TECHNIQUES 

2.1  UTSI  TRANSONIC  WIND  TUNNEL  FACILITY 

Testing  in  conventional  wind  tunnels  usually  results  in  the 
acceptance  of  some  of  the  following  drawbacks:  (1)  scaling 
effect,  i.e.,  limited  Reynolds  number  simulation,  due  (- to  the 
limited  model  size  and  moderate  air  density;  (2)  channel  flow 
effect;  and  (3)  ill-defined  integrated  boundary- layer  thicknesses 
if  the  model  lacks  a clear-cut  leading  edge  or  as  a result  of 
using  a trip  device.  These  drawbacks  will  become  more  significant 
for  flows  at  transonic  speeds. 

In  order  to  avoid  these  drawbacks,  the  present  experiments 
were  conducted  in  a specially  designed  transonic  flow  wind 
tunnel  (see  Refs.  1 and  5).  The  normal  cross-sectional  area  of 
the  test  section  is  12  x 11  inches  (Figure  5) . The  divergence 
angle  between  the  ceiling  and  the  floor  is  approximately  0.8 
degrees  for  the  purpose  of  compensating  for  the  increasing  wall 
boundary  layer  along  the  tunnel  length.  In  order  to  ensure  that 
the  flow  is  free  from  channel  effects  under  different  free-stream 
conditions,  one  adjustable  boundary- layer  suction  device  is  mounted 
on  each  of  the  side  walls  at  94  inches  downstream  from  the  throat 
of  the  nozzle.  These  two  devices  are  also  movable  along  the  flow 
direction  when  necessary.  They  are  designed  to  suck  away  the 
excessive  sidewall  boundary- layer  flow.  The  total  exit  area  of 
the  side  suction  devices  can  be  varied  from  0 to  18  in2.  The 
floor  of  the  wind  tunnel  test  section  was  designed  as  the  test 
model  (see  Ref.  6).  It  has  a total  length  of  144  inches,  with  a 
nominal  width  of  12  inches.  The  leading  edge  of  the  model  is 
located  7 inches  from  the  nozzle  throat.  A natural  boundary 
layer  starting  from  this  leading  edge  is  ensured  by  another 
adjustable  boundary- layer  suction  device  which  is  located  beneath 
the  model  leading  edge.  The  exit  area  of  this  leading-edge 
suction  device  can  be  varied  from  0 to  9 in2.  The  tunnel  floor 
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Figure  5.  Sectional  view  of  The  University  of  Tennessee 
Space  Institute  wind  tunnel. 
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consists  of  a series  of  blocks  which  can  be  modified  or  tailored 
into  any  desired  model  within  the  test  section  length.  The  blocks, 
except  the  7 inch  piece  which  forms  the  leading  edge,  can  be 
arranged  in  any  sequence  and  moved  up  or  down  into  the  test  section 
by  as  much  as  3.75  inches,  as  desired.  The  floor  can  be  used  as 
either  a flat-plate  model,  or  as  a forward- facing  or  a rearward- 
facing step  model  of  different  heights  located  at  different 
stations.  It  can  also  be  arranged  as  a "cavity-like"  model  with 
different  cavity  lengths  and  depths  at  various  locations.  There 
is  enough  space  underneath  the  floor  for  the  installation  of 
various  sensing  instruments.  All  these  wall  pieces  are  movable 
to  facilitate  changes  in  the  model  and  instrumentation  arrange- 
ment. On  top  of  the  test  section  there  is  a 16  inch  high  plenum 
chamber  separated  from  the  test  section  by  a porous  ceiling. 
Ceilings  of  different  porosity  can  be  installed  at  different 
locations  along  the  tunnel  to  minimize  the  wall  interference.  At 
the  end  of  the  test  section  a movable  plug  is  installed  in  the 
diffuser  for  adjusting  the  subsonic  flow  Mach  number  in  the  test 
section. 

2 . 2 MODEL  CONFIGURATIONS 

In  the  experiments  reported  herein,  a forward- facing  step 
model,  five  shallow- cavity  models  and  a flat-plate  model  have 
been  systematically  and  extensively  studied  (Figure  6) . All 
these  models  have  a clear-cut  leading  edge  and  exhibit  two- 
dimensional  flow  features.  The  two -dimensionality  of  flow  was 
tested  by  surface  oil  flow  pattern  observations,  and  by  lateral 
static  and  total  pressure  measurements.  All  of  these  measure- 
ments indicated  that  the  flow  is  in  general  very  close  to  two- 
dimensional  . 

The  forward- facing  step  model  was  formed  by  a 1-inch  high 
normal  step  spanning  the  12-inch  wide  floor  of  the  test  section 
at  74  inches  downstream  from  the  leading  edge  (Figure  6a) . The 
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undisturbed  turbulent  boundary- layer  thickness  was  about  1-inch 
at  this  station.  It  was  tested  to  show  the  Reynolds  number  and 
the  Mach  number  effects  on  flow  separation. 

Shallow- cavity  models  of  equal  separation-  and  recompression- 
step  heights  were  used  to  study  the  phenomenon  of  wake- separation 
interaction  (Figures  6b-f).  These  models  provide  a physical  indi- 
cation of  the  features  of  rear  separation  when  the  oncoming  flow 
was  disturbed  by  a rearward- facing  step.  By  comparing  the  results 
obtained  from  the  "cavity- like”  model  with  that  from  forward- 
facing step  model,  the  effect  of  an  upstream  disturbance  on  the 
downstream  separation  became  clear.  All  cavities  were  maintained 
at  a depth  of  1 inch  and  had  five  different  lengths,  i.e.,  40, 
27.5,  23.5,  15.75,  and  11  inches,  respectively.  These  different 
cavity  lengths  were  used  to  study  the  effect  of  upstream  dis- 
turbance relaxation  and  its  consequence  on  the  downstream  sepa- 
ration. The  question  of  whether  the  separation  streamline  re- 
attaches or  not  for  the  shallow  cavity  could  also  be  examined 
from  the  results  for  different  cavity  lengths.  A flat-plate  model 
was  tested  under  the  same  free-stream  conditions  for  a reference 
of  comparison. 

Surface  pressure  taps,  0.042  inch  diameter,  were  distributed 
longitudinally,  2 inches  to  the  right  of  the  center  line  of  the 
model  surface.  Three  different  spaces  between  taps,  0.25,  0.5  and 
1 inch  respectively,  were  employed  (Figure  7),  These  gave  closer 
measurements  in  regions  where  the  surface  pressure  variation  was 
critical . 

In  order  to  measure  the  total  pressure  variations  across  the 
boundary  layer  (at  any  desired  station  on  the  model)  using  a 
traversing  pitot-probe,  the  floor  block  was  slotted  along  the 
center  line  and  filled  with  small  blocks  of  different  lengths, 
as  shown  in  Figure  8.  The  traversing  pitot-probe  is  screw- 
mounted  such  that  it  can  slide  freely  along  this  lot. 
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b„  40-inch  long  shallow-cavity  model 


c.  27.5-inch  long  shallow- cavity  model 
Figure  6,  Schematic  diagram  of  models. 
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f.  11-inch  long  shallow-cavity  model 
Figure  6.  Concluded „ 
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2 . 3 INSTRUMENTATION 

2.3.1  Traversing  Pitot-Probe 

A variable  speed  vertical- traversing  pitot-probe  was 
designed  to  measure  continuously  the  variation  of  the  total  pres- 
sure across  the  boundary  layer  (Figures  9a  and  b) , The  travers- 
ing mechanism  was  composed  of  a Travel  Unislide-assembly  and  a 
Speed  Control,  with  an  approximate  motor  speed  range  of  1 to  500 
rpm  which  converted  to  a linear  probe  velocity  of  0.0007  to  0.45 
in. /sec.  The  variable  traversing  speed  permitted  a very  slow 
traverse  in  the  inner  viscous  layer  and  an  adequate  speed  in  outer 
region,  so  that  better  measurements  could  be  performed  under  an 
air- saving  running  condition.  A maximum  traversing  distance  of 
3.5  inches  was  possible.  The  probe  was  constructed  from  stain- 
less steel  tubing  with  a squared  end  and  a circular  opening. 

There  is  an  interchangeability  between  three  different  tips  of 
0.04,  0.023,  and  0.011  inches  outside  diameter,  and  0.027,  0.0165, 
and  0.006  inside  diameter,  respectively,  all  of  which  have  been 
used  for  the  measurements  (see  Ref.  1) . However,  the  major  data 
presented  herein  is  based  on  0.023  outside  diameter  probe. 

Since  the  pitot-probe  tip  opening  is  circular  and  its  diameter 
is  much  less  than  4 percent  of  the  boundary-layer  thickness,  the 
data  acquired  is  fairly  reliable,  at  least  in  the  buffer  and  the 
wake  regions,  and  the  displacement  effect  is  also  negligible 
(Refs.  1,  7,  and  8). 

For  the  total  pressure  measurement  in  the  viscous  layer,  the 
flow  turbulence  was  believed  to  be  the  important  source  of 
measuring  error,  especially  in  the  region  very  close  to  the  wall 
®u^fsce.  In  order  to  minimize  the  flow  turbulence  error,  the 
traversing  speed  of  the  pitot-probe  was  limited  to  0.35  in. /sec. 
or  less  through  the  entire  measurements.  Besides  minimizing 
the  flow  turbulence  error,  the  low  traversing  speed  was  important 
in  two  other  aspects.  The  first  is  that  this  speed  was  approxi- 
mately one-millionth  of  the  flow  speed  at  a location  of  0.05 
inches  above  the  model  surface  when  the  local  external  Mach 
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a.  Cross-sectional  view  (schematic) 
Figure  9.  Traversing  pitot  probe  details. 
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Figure 


Photograph 
9.  Concluded. 
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number  was  about  0.6  and  it  was  even  less  than  this  with 
increasing  Mach  number.  Hence,  the  effect  of  the  probe  travers- 
ing speed  on  the  measured  data  was  negligible.  Secondly,  the 
time  needed  to  traverse  a 0.5-inch  thick  boundary  layer  at  this 
speed  was  1.4  seconds.  This  would  render  the  measurements  free 
from  time-lag  effects  (Ref.  6).  Another  difficulty  in  the  total 
pressure  measurement  arose  when  the  measurement  was  taken  close 
to  the  separation  or  the  reattachment  points  where  the  stream- 
lines are  included  to  the  probe  direction.  The  circular  shape 
of  the  tip  employed  was  helpful,  since  the  permissible  yaw  angle 
of  a circular  tip  can  be  as  high  as  ±16  degrees  or  more  with  a 
measuring  accuracy  of  1 percent  (Ref.  8).  At  all  locations 
where  the  pitot  pressure  was  measured,  the  inclined  angle  of  the 
streamlines  was  less  than  16  degrees.  Therefore,  the  error  due 
to  the  misalignment  of  the  pitot- tube  to  the  flow  was  also 
negligible , 

2.3,2  Wall  Static  Pressure  Taps 

The  wall  static  pressures  along  the  model  surface  were 
measured  by  a set  of  static  pressure  taps . The  pressure  taps 
were  formed  by  embedding  a 0,042- inch  inside  diameter  stainless 
steel  tubing  in  each  of  the  0.062-inch  diameter  holes  which  were 
drilled  through  the  model  plate.  One  end  of  the  tubing  was 
flush  with  the  model  surface  and  the  other  end  was  extended  out- 
side the  wind  tunnel  beneath  the  floor  and  connected  to  either 
a manometer  board  or  a pressure  transducer  through  a 0.06-inch 
inside  diameter  flexible  tube. 

The  tap  which  was  located  on  the  model  surface  at  the  same 
distance  from  leading  edge  as  the  tip  of  the  pitot  tube  was 
connected  to  a Statham  pressure  transducer  (PM-131TC;  0 ~ 12.5 
psid) . The  static  pressure  obtained  by  this  tap  was  combined 
with  the  pitot  pressure  to  yield  the  velocity  distribution  across 
the  boundary  layer.  All  other  taps  were  connected  to  a vertical 
manometer  board  which  consisted  of  five  ten-tube  racks,  two  racks 


29 


AE  DC-TR-7  7-1  03 


of  which  were  filled  with  mercury  and  the  others  with  liquid 
of  specific  weight  0.0361  lbf/in3.  The  manometer  permitted  a 
reading  of  0.01-inch  accuracy. 

2.3.3  Skin  Friction  Gauge. 

Some  measurements  of  wind  tunnel  wall  skin  friction  were 
obtained  using  a floating  element  balance  which  was  on  loan 
from  the  Naval  Ordnance  Laboratory,  White  Oak,  Maryland.  The 
balance  is  a null-type  instrument  because  deflections  of  the 
floating  element  due  to  the  wall  shear  force  are  reduced  to  zero 
by  a servomechanism.  The  shear  force  causes  the  balance  arm 
(see  Figure  10)  to  rotate  by  an  amount  proportional  to  the  magni- 
tude of  the  shear  stresses  on  the  surface  of  the  element.  Move- 
ment of  the  balance  arm  is  sensed  by  a translational  Linear 
Variable  Differential  Transformer  (LVDT) . The  LVDT  coil  is 
attached  to  the  balance  housing  while  the  transformer  core  is 
attached  to  the  balance  arm.  Movement  of  the  balance  arm  thus 
produces  a differential  voltage  across  the  output  terminals  of 
the  transformer.  This  voltage  is  converted  into  a plus  or  minus 
DC  voltage  which  is  amplified  and  fed  to  the  servomotor.  Torque 
is  transmitted  from  the  sevomotor  and  across  a gear  train  to  the 
lead  screw.  The  spring  guide  moves  along  the  lead  screw  as  it 
rotates,  producing  a linear  motion  of  the  end  of  the  coil  spring. 
The  force  of  the  spring  acts  to  restore  the  balance  arm  to  a null 
position.  A potentiometer  is  geared  to  the  servomotor  to  provide 
an  analogue  of  the  balance  arm  deflection. 

The  balance  was  mounted  in  the  center  of  the  floor  of  the 
tunnel.  The  floating  element,  of  diameter  0.790  inches,  was 
mounted  in  a hole  of  diameter  0 . 800  inches . The  hole  was  con- 
tained in  a "dummy"  block  which  could  be  moved  to  any  axial 
location  within  the  tunnel. 

The  balance  was  calibrated  by  attaching  weights  to  a fiber 
which  was  fastened  to  the  floating  element  and  suspended  over  a 
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Figure  10.  Balance  schematic. 
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pivot  with  jeweled  bearings.  The  balance  can  measure  forces  up 
to  approximately  10  grams  and  will  null  in  2 to  3 seconds. 

3.0  FLOW  MEASUREMENTS 

3.1  PITOT  PRESSURE  MEASUREMENTS 

Total  pressure  profiles  were  systematically  measured  along 
different  models  at  different  Mach  numbers  and  Reynolds 
numbers.  The  measuring  schedule  is  summarized  in  Table  I. 

Table  I shows  that  the  total  pressure  profiles  along  the 
forward- facing  step  model,  the  40-  and  23.5-inch  shallow-cavity 
models  were  closely  spaced.  The  data  obtained  from  these  models 
will  be  used  as  the  foundation  for  the  study  of  the  wake- 
separation  interaction,  whereas  data  from  the  other  cavities  will 
be  used  to  study  the  effects  of  cavity  length.  Data  from  the 
flat-plate  model  were  used  as  reference  conditions.  The  profile 
at  the  separation  point  was  not  measured  since  the  space  between 
the  separation  point  and  the  forward-facing  step  was  not  long 
enough  to  accommodate  the  probe.  Fortunately,  the  velocity 
profiles  through  the  shear  layer  3 inches  ahead  of  the  separation 
point  were  considered  to  be  sufficiently  downstream  to  show  the 
representative  feature  of  the  flow  close  to  the  step.  These 
profiles  will  be  discussed  later. 

3.2  STATIC  PRESSURE  MEASUREMENTS 

The  static  pressure  distributions  along  the  11- , 15. 75-, 
23.5-,  27.5-,  40-inch  shallow- cavity  models  and  the  forward- 
facing step  model  which  were  used  for  the  wake-separation 
interaction  study  were  obtained  by  using  the  manometer  board. 

The  results  are  shown  in  Figures  11  to  16.  The  pressure  taps 
were  spaced  rather  closely  (0.25-inch  spacing)  in  the  critical 
regions  where  the  pressure  gradients  were  large  or  the  pressure 
gradient  began  to  change  rapidly.  At  locations  where  the  total 
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TABLE  I 

TOTAL  PRESSURE  PROFILE  MEASUREMENT  SCHEDULE 


Model 

Free  Stream 
Mach  Number 

CM  ) 

v oo' 

Unit  Reynolds 
Number 
(Re/ in, ) 

Station 
(inches  from 
leading  edge) 

Flat  Plate 

0.4,  0.5 

0.29  x 

106 

34,  35.5,  49.5 

0.6,  0.7 

to 

51,  61,  64.5 

0.8,  0.9 

0.79  x 

106 

70 

Forward- facing 

0.6 

0.64  x 

106 

20,  47,  51,  61 

Step 

0.  7 

to 

66,  70 

0.75 

0.75  x 

106 

0.9 

40Tr  Cavity 

0.6 

0.44  x 

106 

41,  46,  51, 

61,  66,  70 

27.5"  Cavity 

0.5 

0.39  x 

106 

54,  61,  70 

23.5"  Cavity 

0.8 

0.97  x 

106 

40.5,  41.5, 

42.5,  48.5, 

0.6 

0.44  x 

10  6 

48.5 

15.75"  Cavity 

0.6 

0.44  x 

10* 

44.5 

0.8 

0.97  x 

106 

44.5 

11"  Cavity 

0.6 

0.44  x 

106 

42.5 
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a.  = 0.56,  ReLR  = 14.8  x 1()6,  h - 1 inch 
Figure  11.  Surface  pressure  distribution  along  11-inch  cavity. 


AEDC-TR-77-1  03 


x (in. ) 


0.83,  ReL  = 34.4  x 106,  h 
R 

Figure  11.  Concluded. 


1 inch 


x ( in . ) 


a»  0,56,  ReT  - 1^.9  x 106,  h = 1 inch 

R 

Figure  12.  Surface  pressure  distribution  along  15,75-inch  cavity. 


x (in.) 

0.79,  ReT  = 29.9  x 106,  h 
lR 


= 1 inch 


Figure  12.  Concluded 


00 


^oo  0.58,  Re^  - 15,8  x 106,  h = 1 inch 


Figure  13.  Surface  pressure  distribution  along  23.5-inch  cavity 
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X (in.) 

0.84,  Re^  = 32.9  x 10^,  h = 1 inch 
R 

Figure  13.  Continued. 


x (in. } 


0.85,  ReLR  = 
Figure  13. 


37.0  x 106,  h 
Continued. 


= 1 inch 
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x (in. ) 


0.79,  ReT  = 
Figure  13. 


22.4  x 10^,  h,  = 1 inch 
Continued. 


-p- 


go  M 0,77,  Re,  — 21,5  x 10^,  h=  1 inch 

co  Lj^  ' 

Figure  13.  Continued. 
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0.84, 


x (in.) 

ReT  = 24.8  x 106,  h 
lR 


= 1 inch 


Figure  13.  Continued 
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x ( in. ) 


M =0.60,  ReT  = 15.5  x 106,  h = 1 inch 

“ Lr 

Figure  14.  Surface  pressure  distribution  along  27.5-inch  cavity 


Figure  15.  Surface  pressure  distribution  along  40-inch  cavity. 
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x ( in.  ) 


a.  = 0.61,  pt  - 17.9  psia,  h = 1 inch,  Re^  = 33.1  x 10^ 

Figure  16.  Surface  pressure  distribution  along  forward- facing 
step  model. 
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x ( in . ) 

b»  = 0.91,  pt  = 25  psia,  h = 1 inch,  Re^  = 56.9  x 106 


Figure  16.  Concluded 
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pressure  profiles  were  measured,  the  pressure  taps  were  con- 
nected to  a triple  joint.  The  static  pressure  at  these  locations 
was  sensed  by  a manometer  tube  and  a transducer  simultaneously. 
Data  sensed  by  the  transducer  was  used  to  calculate  the  velocity 
distribution  across  the  shear  layer.  The  difference  between  the 
manometer  and  the  transducer  data  was  less  than  1 percent  as  shown 
in  Figure  17.  The  static  pressure  distribution  on  the  side  wall 
near  the  rearward-facing  step  was  also  measured  (Figure  18) . 

3.3  SURFACE-OIL  FLOW  VISUALIZATION 

Surface  oil  flow  visualization  techniques  have  been  used  to 
extract  details  of  the  separation  and  the  reattachment  regions. 

It  has  also  been  applied  to  the  side  walls  in  order  to  obtain  a 
qualitative  idea  of  the  three-dimensional  structure  of  the  flow. 

A mixture  of  titanium  dioxide  and  kerosene  oil  was  used  for  this 
purpose . 

In  order  to  avoid  an  excessive  oil  accumulation  near  the 
separation  line  and  to  keep  the  oil  from  forming  U-shaped  waves, 
different  proportions  of  titanium  dioxide  to  kerosene  oil  were 
tested  and  used  for  the  different  skin  friction  regions  (Ref.  7) . 
On  the  side-wall,  the  viscosity  of  the  applied  mixture  was  high 
enough  so  that  the  patterns  were  believed  to  be  unaffected  by 
the  gravitational  force. 

A representative  surface  oil  flow  photograph  for  the  flow 

near  the  forward-facing  step  in  the  27.5-inch  cavity  (at  M - 0.6 

6 00 
and  ReL^  = 32.5  x 10  ) is  shown  in  Figure  19.  Data  of  the 

separation  and  reattachment  distances  on  different  models  at 
different  Mach  numbers  and  Reynolds  numbers  obtained  from  the  oil 
flow  observations  are  shown  in  Table  II  and  Figure  20. 

Based  on  this  data,  several  features  of  the  flow  can  be  dis- 
cussed as  follows.  First,  for  the  forward -facing  step  at  the 
end  of  the  cavity: 


51 


Surface  pressure  sensed  by  manometer  tube 

(psia) . 

Figure  17 „ Comparison  of  pressures  sensed  by  pressure 
transducer  and  by  manometer  tube. 


Figure  19.  Photograph  showing  oil  flow  pattern  ahead  of  the  forward- 
facing step  at  = 0.6  and  ReLf  = 32.5  x 106. 
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TABLE  II 

REATTACHMENT  AND  SEPARATION  DISTANCES  FOR  FLOW  OVER  CAVITIES 


Cavity 

Length 

(in.) 

M 

Pt 

(psi.a) 

Reattachment 
Distance,  £ 
(in.)  r 

Separation 
Distance,  £ 
tin,)  s 

% 

V ■*  ■*  > 

11.0 

0.83 

28.2 

7.0 

1.0 

28,44  x 106 

11.0 

0,56 

18,0 

7,0 

1.0 

14.33  x 106 

15.75 

0,80 

25.2 

6,5 

0.9 

24.99  x 10 6 

15,75 

0,66 

25.7 

6,5 

1.0 

22. 94  x 106 

23.5 

0.85 

29.6 

6,0 

1.2 

39,19  x 106 

23.5 

0,84 

26,6 

6,5 

0.9 

27,09  x 106 

23.5 

0.58 

18.6 

5,5 

0.9 

15.26  x 106 

27.5 

0,70 

17.6 

5.2 

1,1 

22,10  x 106 

27.5 

0.60 

17.5 

5.2 

1,1 

20,00  x 106 

27.5 

0.60 

17.7 

5.2 

1.1 

20,23  x 106 

40.0 

m 

17,8 

5.9 

1.1 

17.65  x 106 

40.0 

B9 

16.9 

5.9 

1.1 

15.37  x 106 

40.0 

m 

17.5 

5,9 

1.1 

15,77  x 106 
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Figure  20 


Effect  of  cavity  length  on  reattachment  and 
separation  distance  (M  = 0.56  to  0.85, 

Re/ft  = 5.04  x 106  to  13.8  x 10&) . 
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1.  It  was  seen  that  the  three-dimensional  effects  at  the 
corner  of  the  step  and  the  tunnel  wall  made  the  flow 
pattern  quite  complicated  (Figure  19) . However,  the 
spanwise  influence  of  this  effect  was  not  sufficient 
to  destroy  the  two-dimensional  nature  of  the  flow  over 
the  major  part  of  the  span. 

2.  These  secondary  flows  were  confined  within  a region 
less  than  two  step  heights  from  the  side  walls. 

3.  There  was  a vortex  in  the  flow  ahead  of  the  step  with 
its  center  about  1.8  inches  upstream  of  the  step  and 
0.25  inches  from  the  side  wall. 

4.  The  oil  flow  on  top  of  the  step  indicated  a reattach- 
ment line  of  the  separated  flow  over  the  step  at  about 
2.4  inches  downstream  of  the  step.  There  was  also  an 
accumulation  of  oil  right  at  the  step  comer.  This 
would  result  from  the  reversed  flow  in  the  separation 
bubble  on  the  top  of  the  step  forcing  the  oil  to  the 
foremost  portion  of  the  vortex. 

5 . The  separation  line  on  the  floor  ahead  of  the  forward- 
facing step  was  clearly  seen.  The  distance  from  the 
separation  point  to  the  root  of  the  step  (£)  was  roughly 
constant  (about  one  step  height)  in  the  Mach  number 
range  of  0.6  to  0.8.  It  decreased  slightly  with  in- 
creasing cavity  length  (JL)  when  the  cavity  length  was 
long  enough  and  reached  a minimum  (about  0.9  step 
heights)  at  a cavity  length  of  about  24  inches,  then 
increased  slightly  with  a further  decrease  in  cavity 
length. 

6.  For  the  rearward-facing  step  at  the  beginning  of  the 
cavity  the  following  features  appeared:  the  reattach- 
ment region  behind  the  rearward- facing  step  was  roughly 
defined  by  the  oil  flow.  Several  points  were  clear: 
the  reattachment  flow  was  quite  three-dimensional;  the 
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effect  of  the  cavity  length  on  the  reattachment  dis- 
tance (&r)  had  the  same  trend  as  but  slightly  stronger 
than  that  on  the  separation  distance  ahead  of  a for- 
ward facing  step;  the  minimum  reattachment  distance 
was  about  5.2  step-heights  occurring  at  Mm  ~ 0.6  with  a 

cavity  length  of  27.5  inches  and  the  maximum  was  about 

7.0  step-heights  at  . 0.6  and  0.8  with  a cavity 
length  of  11  inches  (Table  II,  page  45).  It  appears 
that  the  reattachment  distance  is  a function  of  the 
cavity  length  and  almost  independent  of  the  subsonic 
free- stream  Mach  number. 

7.  Within  the  recirculating  region  behind  the  rearward- 
facing step,  there  was  an  oil  line  about  1 inch  from 
the  step  which  indicated  a separation  of  the  reversed 
flow  and  the  presence  of  a secondary  vortex. 

4.0  DISCUSSION  OF  EXPERIMENTAL  RESULTS 

4.1  SURFACE  PRESSURE  DISTRIBUTIONS 

4.1.1  Dimensional  Analysis  and  Reference  Parameters 

In  these  experiments,  it  was  assumed  that  no  heat  transfer 
or  chemical  reactions  took  place.  The  non-dimensional  para- 
meters that  may  be  pertinent  to  the  analysis  of  the  experimental 
data  can  be  obtained  by  applying  a dimensional  analysis.  With 
the  understanding  that, 

p = p (flow  variables  and  model  configurations) 

“ P (Pref'  u’  p'  a'  xi’  gi^ 
where 

Pref  = reference  pressure,  u = flow  velocity,  a = sound  speed, 
and 

x*  = streamwise  characteristic  length,  such  as  cavity  length,  etc. 
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y±  = characteristic  length  perpendicular  to  the  model  surface 
such  as  boundary- layer  integral  thickness,  etc. 


It  results  in: 


pux . y . \ 

1 u J \ 


ref 


a x 


it 


or 


/ 


Pref 


= f 


pux . 


u 

a1 


yi 


0 


(4.1) 


f Re,  M, 


(4.2) 


This  relationship  indicates  the  dominating  parameters  involved 
in  this  study. 

The  significant  static  pressures  concerned  in  this  study 
could  be  the  base  pressure  (pb) , the  pressure  at  the  reattach- 
ment point  (pf) , the  plateau  pressure  (pp) , the  pressure  at  the 

separation  point  (p  ) or  the  peak  pressure  (p  ) (see  Figure  4). 

s tn 

The  most  generally  used  reference  pressures  are  the  total  pres- 
sure (pt)  and  the  free-stream  static  pressure  (p^) . In  order  to 
understand  the  trend  of  the  energy  consumption  by  the  turbulent 
mixing  and  the  recovery  of  the  flow  along  the  model  surface, 
some  other  pressures,  such  as  the  base  pressure  and  the  plateau 
pressure,  are  also  worth  using  as  references. 

The  characteristic  length  (xA)  considered  will  be  the  dis- 
tance from  the  leading  edge  along  the  strearawise  direction  (x, 

Lr  or  Lf) , which  indicates  the  historical  effect  of  the  flow 
turbulence  structure.  Another  streamwise  characteristic  length 
will  be  the  cavity  length  (£)  which  can  partially  indicate  the 
strength  of  the  mutual  interference  effect  between  the  upstream 
disturbance  and  the  downstream  separation.  The  characteristic 
length  perpendicular  to  the  model  surface  (yi)  will  be  the  step 
height  (h)  and  the  boundary- layer  integral  thicknesses  (6,  <5* 
and  6).  The  first  characteristic  y ^ , namely  h,  will  relate  to 
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the  strength  of  the  upstream  disturbance  as  well  as  to  the  down- 
stream separation.  The  last  three  will  indicate  the  state  of  the 
flow  energy  in  the  shear  layer.  These  are  all  very  useful  in 
analyzing  the  features  of  the  flow. 

4.1.2  Surface  Pressure  Distribution  Over  the  Entire  Model 

A general  description  and  discussion  of  the  pressure  distri- 
bution over  the  entire  cavity  model  is  made  here  (see  Figure  21) . 
The  detail  features  of  the  important  pressures  such  as  the  base 
pressure  pfa , the  reattachment  pressure  pr,  the  cavity  plateau 
pressure  pp,  the  separation  pressure  pg,  and  the  peak  pressure 
Pm  and  the  estimation  of  the  critical  cavity  length  (see  p.  65 
for  definition)  based  on  the  features  of  the  pressure  distribution 
over  the  shallow-cavity  model,  as  well  as  mutual  relationships 
among  the  important  pressures  will  be  discussed  separately  in  the 
latter  subsections. 

4. 1.2.1  Flow  Over  Cavities  of  Different  Length  to  Depth  Ratio 

The  surface  pressure  distributions  over  the  cavity-like 
models  of  length  to  depth  ratios  £/h  = 11,  15.75,  23.5,  27.5 
and  40  at  Mro  = 0,6  and  Re^  = 15  x 106  (based  on  the  length 

measured  from  the  leading  edge  to  the  rearward-facing  step  station, 
Lr)  has  been  shown  earlier  in  Figures  11  - 15,  as  p/pt  versus 
x (measured  in  inches  from  the  leading  edge),  respectively. 

These  distributions  are  combined  together  as  the  pressure  co- 
efficient  (C  = (p  - p )/  (%yp  M ) versus  the  non-dimensional  dis- 
tance,  xR/ J,  in  Figure  21.  They  give  a general  picture  of  the 
cavity  length  to  depth  ratio  effect  at  subsonic  speeds. 

From  Figure  21  it  is  found  that  the  pressure  coefficient 
immediately  downstream  of  the  rearward- facing  step  shows  a 
significant  drop  which  increases  as  the  cavity  is  shortened. 

This  is  a typical  behavior  of  the  base  pressure  problem.  The 
pressure  coefficient  drop  nearly  disappers  at  Jl/h  = 11.  The 
location  of  the  minimum  pressure  coefficient  moves  as  the  cavity 
length  or  pressure  coefficient  drop  decreases.  The  pressure 
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coefficient  drop  is  followed  by  a rapid  rise  which  is  steepest 
for  the  longer  cavity  models.  The  mutual  interaction  of  the 
rearward  and  forward  facing  steps  is  very  strong  as  evidenced  by 
the  great  change  in  the  base  pressure  or  the  separation  pressure 
even  though  both  steps  are^up  to  30  to  40  boundary  layer  thick- 
ness distance  apart. 

For  the  longer  cavities  (&/h  > 15.75),  there  is  a pressure 
plateau  along  the  center  portion  of  the  floor,  but  the  plateau 
disappears  as  the  cavity  shortens,  as  in  the  case  of  Jl/h  = 11. 

The  starting  point  of  the  plateau  moves  downstream  as  the  cavity 
length  is  decreased  and  eventually  coincides  with  the  pressure 
peak.  In  the  middle  portion  of  the  cavities,  the  pressure  co- 
efficient increases  at  first  with  shortening  cavity,  but  eventu- 
ally begins  to  decrease.  In  the  center  portion  of  the  cavity, 
the  pressure  coefficient  reaches  a maximum  value  of  about  0.3 
(corresponding  to  £/h  = 15.75),  whereas  the  maximum  peak  pressure 
coefficient  ahead  of  the  forward- facing  step  is  about  0.54 
(corresponding  to  fe/h  = 23.5).  These  values  are  about  50  and  20 
percent  higher  than  those  for  the  comparable  incompressible  flow 
(Ref.  10).  Upstream  of  the  rearward- facing  step,  i.e.,  flow 
ahead  of  the  base,  the  pressure  coefficient  shows  a sinous 
feature.  Downstream  of  the  forward-facing  step  there  is  another 
pressure  coefficient  drop  which  is  much  stronger  than  that  of  a 
flow  over  a rearward-facing  step  and  increases  with  decreasing 
cavity  length.  This  drop  recovers  rapidly  towards  the  free- 
stream  condition. 

The  pressure  distribution  at  transonic  speeds  along  the 

cavity-like  models  of  length  to  depth  ratios  £/h  = 15.75  and 

23.5  at  = 0. 8 and  Re^  - 29  x 10^,  shown  in  Figures  12b  and 

R 

13c  as  p/pt  vs.  x,  are  replotted  in  Figure  22  as  C versus 
xr/£.  The  general  trend  of  the  pressure  coefficient  distributions 
for  these  lower  transonic  speed  flows  is  similar  to  that  of  the 
subsonic  case  shown  in  Figure  21.  But  there  are  several  special 
features  in  the  transonic  flow  cases.  These  are  as  follows: 
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First,  the  pressure  drop  at  the  base  of  these  two  length  to 
depth  ratio  cavities  is  nearly  the  same.  Secondly,  the  peak 
pressure-coefficient  difference  between  models  of  1/ h = 23.5 
and  15.75  is  greater  than  that  of  the  corresponding  subsonic 
case.  Thirdly,  variation  of  the  pressure  coefficient  distri- 
btion  upstream  of  the  rearward- facing  step  is  stronger  than 
that  in  subsonic  flow. 

It  can  be  pointed  out  further  by  observation  of  Figure  23 
that  the  pressure  coefficients  along  the  center  portion  of 
model  of  If h = 23.5  are  coincident  for  the  higher  (M  =0.80 
and  Re^  = 28.6  x 10 6)  and  the  lower  (Mu  = 0.58  and  Re, 

R g 

15.8  x 10  ).  The  coincidence  seems  true  for  all  flows  over 
the  model  of  2,/h  = 23.5  at  different  Mach  numbers  and  Reynolds 
numbers,  which  can  be  verified  by  the  observation  of  Figures 
23,  24,  and  25.  On  the  other  hand,  for  a shorter  model 
(S7h  » 15.75),  the  pressure-coefficient  gradient  along  the 
center  portion  of  a cavity  for  a higher  inertia  energy  flow 
(higher  Reynolds  number)  is  much  lower  than  that  for  a lower 
inertia  energy  flow  and  almost  coincides  with  the  pressure- 
coefficient  gradient  along  the  model  of  Jl/h  = 11.00  at  the 
lower  Reynolds  number  (Figure  26).  The  pressure  plateau  of  the 
higher  inertia  energy  flow  disappears  at  £/h  = 15.75  instead  of 
at  £/h  = 11  as  in  the  low  inertia  energy  case  (Figures  26  and 
21) • It  seems  that  the  effect  of  inertia  energy  upon  the  pres- 
sure coefficient  distribution  becomes  more  sensitive  when  Jl/h  is 
smaller.  And,  it  can  be  clarified  further  in  the  next  subsection 
that  the  inertia  energy,  which  will  affect  the  pressure  coef- 
ficient distribution,  is  mainly  due  to  the  Reynold  number  change. 
For  longer  cavities  (£/h  ^ 23.5)  the  values  of  the  nearly  con- 
stant pressure  over  the  midsection  of  the  cavity  are  approximately 
the  mean  value  of  the  corresponding  maximum  and  minimum  pressures 
in  the  cavity.  In  other  words,  the  magnitude  of  the  pressure 
drop  over  the  rearward- facing  step  is  approximately  equal  to  the 
pressure  rise  induced  by  a forward- facing  step  of  the  same  height. 
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Figure  24.  Effect  of  Reynolds  number 
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Figure  25.  Effect  of  Mach  number 
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However,  as  if h becomes  small  (fc/h  ^ 15.75),  the  pressure  rise 
ahead  of  the  forward-facing  step  becomes  less  than  the  pressure 
drop  over  the  back-step. 

From  the  above  description,  it  is  found  that  the  general 
trend  of  the  surface  pressure  distribution  over  the  cavity 
floor  of  different  depth  to  length  ratios  is  similar  to  the  incom- 
pressible flow  case.  However,  it  is  not  possible  to  predict  the 
pressure  coefficient  of  a compressible  flow  over  a cavity-like 
model  from  a corresponding  incompressible  result  by  using  any 
existing  rules.  This  is  the  case  even  in  the  region  where  the 
separation  and  reattachment  regions  are  excluded.  For  instance, 
consider  using  the  well-known  Prandtl-Glauert  rule  to  simulate 
the  flow  at  Mot  = 0.79  from  the  measured  flow  at  = 0,19.  The 
result  shows  a tremendous  error  when  compared  with  the  measured 
data.  In  the  center  portion  of  the  cavity,  the  error  of  the 
corrected  pressure  coefficient  can  be  as  large  as  50  percent 
(Figure  27) . 

4. 1.2.2  Flow  Over  Cavity  With  Different  Reynolds  Numbers  and 
Mach  Numbers . 

Figure  24  shows  the  pressure  coefficient  distribution  for 
flow  over  a 23.5  inch  cavity-like  model  at  - 0.8  and  different 
Reynolds  numbers  (22.4  x 106  and  28,6  x 106),  It  can  be  seen 
from  this  plot  that  at  this  transonic  speed,  the  higher  Reynolds 
number  flow  will  cause  a lower  base  pressure  and  a higher  peak 
pressure,  whereas  the  plateau  pressure  is  relatively  insensitive 
to  the  Reynolds  number  change.  Figure  25  shows  that  the  pres- 
sure coefficient  distribution  seems  insensitive  to  the  variation 
of  Mach  number  along  the  entire  model. 

Figures  23  and  24  give  evidence  that  the  flow  over  a 
cavity-like  model  will  have  a lower  base  pressure  and  a higher 
peak  pressure  if  the  inertia  energy  of  the  flow  is  higher. 

Being  aware  of  the  facts  described  above,  it  can  be  further 
concluded  that  this  phenomenon  is  mainly  attributed  to  the 
Reynolds  number. 
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Figure  27.  Comparison  of  measured  data  and 
Prandtl-Glauert  rule. 


AEDC-TR-77-103 


AEDC-TR -77-103 


4. 1.2.3  Pressure  Distribution  Over  Forward-Facing  Steps, 

The  pressure  distributions  over  a forward-facing  step  model 
are  given  in  Figures  16a  and  b,  as  p/pc  versus  x,  and,  in 
Figure  28  as  p/p^  versus  xf/h,  A comparison  of  the  pressure 
distribution  for  turbulent  flow  over  a forward- facing  step  model 
at  supersonic,  transonic,  and  subsonic  speeds  is  shown  in  Figure 

29.  These  figures  show  that  in  the  lower  transonic  and  subsonic 
ranges,  the  pressure  distributions  over  a forward-facing  step 
model  have  a peak  pressure  point  about  one  step  height  upstream 
of  the  step,  no  matter  how  different  the  Mach  numbers  and 
Reynolds  numbers  are.  The  upper  transonic  and  supersonic  flows 
have  a pressure  plateau  in  addition  to,  and  ahead  of,  the  pres- 
sure peak,  and  the  supersonic  flow  reaches  the  highest  peak 
pressure  and  has  the  shortest  upstream  influence  distance. 

In  subsonic  flows,  the  pressure  rise  ahead  of  the  step  is  simply 
due  to  the  conversion  of  the  kinetic  energy  of  the  flow  into 
pressure  energy  by  the  presence  of  the  step.  It  is  natural 
that  the  flow  having  the  higher  speed  will  have  a higher  pres- 
sure peak  (if  the  other  conditions  remain  the  same) . As  the 
flow  speed  becomes  greater  than  sonic,  a shock  wave  is  formed 
due  to  the  turning  of  the  dividing  streamline  ahead  of  the  step. 
This  shock  and  the  inviscid  flow  field  result  in  the  so-called 
plateau  pressure  regime.  As  the  step  is  approached,  it  forms 
the  peak  pressure  because  the  flow  is  required  to  slow  down.  By 
observing  the  difference  in  subsonic  and  supersonic  flow  sepa- 
ration, it  seems  that  the  original  separation  pressure  transfers 
to  the  plateau  pressure  while  the  separation  point  suddenly 
moves  in  the  upstream  direction  as  the  approaching  Mach  number 
increases  from  below  to  over  one.  This  fact  is  shown  in  Figure 

30.  In  other  words,  if  the  free-stream  Mach  number  decreased 
from  supersonic  to  subsonic  flow,  the  distinguished  plateau  and 
separation  pressures  merge  together  to  become  a subsonic  flow 
separation  pressure  (Figure  30) . 
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Figure  30,  Effect  of  free-stream  Mach  number  on  the 
separation  and  plateau  pressures. 
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4.1.2. 4 Discussion  of  Critical  Cavity  Length  Eased  on  Pressure 
Measurements 

When  the  fluid  flows  over  a cavity,  the  flow  will  separate 
from  the  corner  of  the  rearward-facing  step  first,  and  then  re- 
attach to  the  floor  of  the  cavity  if  the  length  to  depth  ratio 
of  the  cavity  is  large  enough.  The  flow  then  reseparates  as  it 
approaches  the  forward-facing  step.  This  type  of  flow  is  called 
flow  over  a "closed"  cavity.  It  is  commonly  known  that  two 
major  circulating  regions  exist  in  this  type  of  flow,  in  the  re- 
attaching bubble  at  the  rearward-facing  step  side  and  in  the  re- 
separation bubble  at  the  downstream  forward-facing  step  side, 
respectively.  The  distance  between  the  reattachment  point  and 
the  downstream  separation  point  in  the  cavity  decreases  as  £/ h 
is  decreased.  As  the  length  to  depth  ratio  decreases  to  a certain 
value,  the  flow  becomes  unsteady--the  flow  will  reattach  to  the 
cavity  floor  or  jump  from  the  rearward-facing  step  to  the  down- 
stream forward- facing  step  alternately.  Further  decrease  of  J4/h 
will  cause  the  flow  to  jump  firmly  across  the  rearward-facing 
step  to  the  downstream  forward- facing  step.  In  such  a case  no 
reattachment  point  can  be  found  on  the  floor,  and  the  flow  be- 
comes steady  again.  This  type  of  flow  is  called  flow  over  an 
"open"  cavity  because  the  dividing  streamline  has  jumped  across 
the  cavity  and  never  reattaches  to  the  floor. 

The  pressure  distribution  in  a cavity  shows  some  system- 
atical change  corresponding  to  the  variation  of  the  length  to 
depth  ratio  (Figure  21) . This  change  can  be  characterized  by  the 
pressure  measurement  at  the  midpoint  on  the  floor  of  the  cavity 
(Cp£/2) . Figure  31  shows  that  the  pressure  at  the  midpoint  of  the 
floor  increases  with  the  initial  h/£  increase,  It  reaches  a 
maximum  value  and  then  decreases.  This  maximum  value  corresponds 
to  the  envelope  value  of  the  Cp  distribution  in  the  center  portion 
of  cavities  with  different  length  to  depth  ratios  (Figure  21) . 

It  is  found  in  Figure  31  that  (C  „ .„)  for  a subsonic  flow  is 

° p £ / 2 max 
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Figure  31.  Effect  of  depth  to  length  ratio  on  the 
pressure  coefficients  at  the  midpoint 
of  a cavity. 
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higher  than  that  of  the  corresponding  incompressible  flow  case. 

The  only  C data  available  at  the  supersonic  speed  shows  a 
lower  value  than  that  for  the  incompressible  flow.  Hence  there 
is  no  consistent  increase  in  this  pressure  as  the  free-stream 
velocity  is  increased  through  subsonic  values  into  the  super- 
sonic range.  The  change  of  slope  of  C , in  Figure  31  from 

P X / 4- 

positive  to  negative  with  increasing  h/ a must  correspond  to  some 
drastic  change  of  the  flow  structure  in  the  cavity.  This  point 
of  view  can  be  verified  by  the  total  pressure  measurements 
(Cpt  = (pt-P£/2)/(%TPA/2MJ>  with  a pitot-probe  placed  on  the 
floor  at  the  midpoint  of  the  cavity.  The  results  show  that  Cp)_ 

decreases  with  increasing  h/ i and  reaches  a zero  value  at  h /i 
0.095  (i  = 10.5  in.)  which  is  slightly  larger  than  the  value 
corresponding  to  the  maximum  C , The  zero  value  of  C re- 

p£/  2 p 

veals  that  a separation  occurs  at  the  midpoint  of  the  floor,  which 
is  evidence  of  an  open  cavity.  Therefore,  it  can  be  concluded  that 
the  opening  of  the  cavity  will  be  impending  when  the  cavity  is 
being  shortened  continuously  until  the  envelope  value  of  the 
pressure  distribution  in  the  center  portion  of  the  cavity  is 
reached. 

Figure  31  shows  that  the  subsonic  flow  (M  = 0.6,  Re  = 

6 “ R 
15.9  x 10  ) has  a shorter  critical  cavity  length  that  the  incom- 
pressible flow  has.  Figure  26  shows  that  a transonic  flow  ' 

= 0.79,  Re,  = 29.9  x 106)  over  a cavity  of  £/h  = 15.75  has 
R 

the  same  pressure  distribution  along  the  center  portion  of  the 
cavity  as  a subsonic  flow  (M  =0.56,  ReT  = 15.8  x 10 6)  over  a 

CO  ij  ' 

R. 

cavity  of  £/h  = 11.  This  phenomenon  indicates  that  at  this  tran- 
sonic speed  the  critical  cavity  length  might  be  longer  than  that 
of  the  subsonic  flow  case.  But,  as  Mach  number  increases  further 

(M  =0.83,  Re,  = 34.4  x 10^),  the  critical  cavity  length  de- 
R 

creases  again  (Figure  32) . It  seems  that  the  critical  cavity 
length  will  be  a maximum  at  some  particular  transonic  speed  and 
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Figure  32 „ Effect  of  flow  Inertia  energy  on  the  pressure 
coefficient  distribution  over  a shallow  cavity 
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further  tests  are  needed  to  clarify  this  point.  The  flow  might 
be  very  non-steady  when  the  cavity  flow  is  transferring  from  the 
closed  to  the  open  mode.  Therefore,  it  is  actually  impossible 
to  define  a single  value  for  the  critical  cavity  length  corres- 
ponding to  a given  flow  condition.  Instead,  a range  of  critical 
cavity  lengths  would  be  estimated  for  different  Mach  numbers  as 
shown  schematically  in  Figure  33. 

4.1.3  Discussion  of  Base  Pressure 

It  can  be  found  in  Figure  34  that  for  flow  over  a rearward- 
facing step  at  high  subsonic  of  transonic  speeds  with  a separation 
bubble  at  the  downstream  side,  the  non-dimensional  base  pressure 
(Pb/Pj  decreases  linearly  with  an  increase  in  the  Reynolds 

number  up  to  Re  = 30  x 106,  At  higher  Reynolds  numbers,  the 

R 

Pb/Pro  versus  Re^  curve  keeps  the  same  slope  but  shifts  to  a 

R 

relatively  higher  value.  Since  there  are  only  two  data  points 
(M^  = 0.84  and  0.85  respectively)  in  the  region  of  ReT  > 30  x 106 

LR 

and  only  one  data  point  (M  = 0.84)  in  the  region  of  ReT  < 30  x 
10  , what  will  happen  for  flow  speeds  of  M x 0.8  and  M >.0.84 

00  oo  ' 

in  the  range  of  Re^  > 30  x 106  and  Re.  < 30  x 106  , respectively, 

R lr 

needs  to  be  clarified  by  further  tests. 

The  same  data  which  are  shown  in  Figure  34  are  replotted  as 

0^  versus  M^  in  Figure  35.  For  flow  at  a Mach  number  less  than 
b 

0.8,  the  pressure  coefficient  distribution  increases  slightly  with 
Mach  number  at  a constant  total  pressure.  As  the  Mach  number 
becomes  greater  than  0.8,  the  pressure  coefficient  changes  very 
peculiarly . 

Figures  36  and  37  show  the  influence  of  the  downstream 

separation  on  the  upstream  disturbance.  The  non-dimensional  base 

pressure  (p  /p  or  C ) increases  with  h/i,  This  means  that  a 
U . Pb 

separation  bubble  downstream  with  a certain  given  strength 
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Figure  33.  Nondimensional  critical  cavity  length  (schematic 
representation  of  possible  result). 


Figure  34.  Effect  of  Reynolds  number  on  the  base  to  free- 
stream  pressure  ratio. 
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Figure  35 „ Effect  of  free-stream  Mach  number  on  the  base 
pressure  coefficient. 
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Figure  36,  Effect  of  depth  to  length  ratio  on  the  base 
to  free-stream  pressure  ratio. 
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Figure  37,  Effect  of  depth  to  length  ratio  on  the  base 
pressure  coefficient. 
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(Moo>  Re,  , Z and  h)  will  enhance  the  base  pressure  indicated  at 
. R 

an  upstream  location  when  the  separation  bubhle  moves  towards  the 

base.  As  h/ Z approaches  zero,  C will  reach  a negative  con- 

ph. 

stant  or  zero  depending  upon  whether  Z+o  or  h+o,  Whether  or  not 
the  C verus  h/ z curve  for  flow  over  cavities  of  different 

*D 

lengths  and  a fixed  step  height  will  coincide  with  that  of  dif- 
ferent step  heights  and  a fixed  cavity  length  needs  to  be 
clarified  by  further  tests. 

In  the  discussion  of  the  base  pressure,  the  approaching 
boundary- layer  thickness  (6,  ) at  the  top  corner  of  the  rearward- 
facing  step  is  usually  a significant  parameter.  It  represents 
the  combined  effect  of  the  free-stream  Mach  number  CM  ) . the 
Reynolds  number  (Re^  ) and  the  historical  memory  of  the  surface 

roghness  and  other  aspects  of  the  upstream  flow  condition.  In 

the  present  study,  may  also  show  the  effect  of  different 

R 

downstream  adverse  pressure  gradients  which  are  more  or  less 
severe  depending  on  the  depth  and  the  length  of  the  cavity. 

Therefore,  the  measurement  of  $ and  a plot  of  p versus  5 may 

R b Lr 

be  worthwhile  in  a further  study. 

4.1.4  Discussion  of  Reattachment  Pressure 

In  the  reattachment  zone,  the  pressure  recovers  almost 
immediately  in  comparison  with  the  long  relaxation  process  of 
the  turbulence  and  velocity  profiles  (see  Section  4.3).  The 
conditions  at  the  reattachment  point  are  controlled  by  a mass 
balance  in  the  so-called  dead-air  region  (but  the  velocities  in 
this  region  can  be  quite  high).  To  accommodate  this  process, 
the  pressure  rise  to  reattachment  must  be  sufficient  to  return 
that  quantity  of  mass  entrained  by  the  shear  layer  back  to  the 
dead-air  region.  The  reattachment  point  is  the  point  where  the 
flow  bifurcates  and  a new  viscous  sublayer  starts  downstream. 
Properties  at  certain  points  sufficiently  far  downstream  of  the 
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reattachment  point  are  important  for  rear  separation  studies 

since  all  historical  effects  upon  the  boundary  layer  hitherto 

are  embedded  in  the  flow  properties  at  this  point.  These 

properties  form  a new  upstream  boundary- layer  condition  for  the 

corresponding  rear  separation.  The  reattachment  pressure  is  the 

static  pressure  at  the  downstream  end  of  the  dividing  streamline. 

The  compression  process  will  recover  the  base  pressure  to  a 

plateau  value  higher  than  the  reattachment  pressure  and  equal  to 

or  even  higher  than  the  corresponding  free-stream  pressure.  At 

supersonic  speeds,  the  magnitude  of  the  reattachment  pressure  is 

roughly  the  mean  value  of  the  corresponding  base  and  plateau 

pressues,  whereas  at  subsonic  or  transonic  speeds,  it  is  about 

2/3  of  the  base  pressure  value.  Also,  as  it  has  already  been 

pointed  out,  for  supersonic  speeds  the  reattachment  distance  is 

much  shorter  than  for  the  subsonic  case.  From  these  two  pieces 

of  evidence,  it  can  be  concluded  that  the  supersonic  flow  has  a 

greater  ability  for  pressure  recovery.  By  observation  of  Figure 

38,  it  is  found  that  the  reattachment  pressure  increases  linearly 

with  the  Reynolds  number  (ReT  ) at  a given  liach  number,  while  at 

lR 

a given  Reynolds  number,  the  reattachment  pressure  decreases 
linearly  with  increasing  Mach  number  (Figure  39).  Figures  40  and 
41  show  that  C and  Pr/poo  increase  with  h/£.  It  is  also  found 

in  Figure  40  that  for  h/ Z = 0.025  (Jl  = 40  in.)  the  magnitude  of 
reattachment  pressure  coefficient  is  less  than  zero.  Theoretically, 
there  are  two  possibilities  as  h/£  approaches  zero:  first,  if 
the  cavity  length  is  kept  constant,  the  reattachment  pressure  co- 
efficient will  be  zero  as  the  cavity  depth  approaches  zero. 

Second,  if  cavity  depth  is  a constant,  the  reattachment  pressure 
coefficient  will  approach  the  value  corresponding  to  the  flow 
over  a rearward-facing  step  without  a rear  separation,  where  the 
pressure  coefficient  is  usually  less  than  zero,  as  the  cavity 
length  exceeds  a certain  value.  Data  in  Figures  40  and  41  were 
obtained  by  fixing  the  step  height  at  h = 1 inch  and  varying  the 
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Figure  38.  Effect  of  Reynolds  number  on  the  reattachraent 
to  free-stream  pressure  ratio. 
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Figure  39.  Effect  of  free-stream  Mach  number  on  the 
reattachment  pressure  coefficient. 
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cavity  length  from  11  inches  to  4Q  inches.  If  the  test  is  con- 
ducted by  fixing  the  cavity  length  and  varying  the  step  height, 
different  results  may  be  obtained,  at  least  in  the  range  for 
which  h/ l approaches  zero, 

4.1.5  Discussion  of  the  Cavity  Plateau  Pressure 

Cavity  plateau  pressure  is  also  a function  of  the  Reynolds 
number  but  the  dependency  is  weak  as  compared  with  the  base 
pressure-Reynolds  number  relation  (Figures  34  and  42),  The 
trend  of  dependence  of  the  non-dimensional  plateau  pressure 
(Pp/pJ  and  the  base  pressure  (pfa /p^)  on  Reynolds  number  are 
in  the  opposite  sense.  The  former  shows  a slight  increasing 
trend  with  Reynolds  number  whereas  the  latter  decreases  with  in- 
creasing Reynolds  number.  The  effect  of  Mach  number  upon  the 
non-dimensional  plateau  pressure  p / p is  to  increase  p /p  with 

p*cn  ^ 5 ' r cy 

increasing  Mach  number  (Figure  43).  Flow  over  a rearward- facing 
step  can  be  recompressed  from  the  base  pressure  to  a higher 
plateau  value  by  increasing  either  the  Reynolds  number  (Figure  42) 
or  the  free-stream  Mach  number  (Figure  43)  or  by  introducing  an 
adverse  pressure  gradient  downstream  (Figure  44) . Since  the 
magnitude  of  the  cavity  plateau  pressure  is  larger  than  the  free- 
stream  pressure,  the  rearward- facing  step  performs  somewhat  like 
a vortex  generator  and  acts  to  suppress  the  downstream  separation. 
The  efficiency  of  the  energizing  process  will  increase  with 
Reynolds  number,  free-stream  Mach  number  or  with  the  depth  to 
length  ratio  of  the  cavity.  From  the  wake-separation  interaction 
point  of  view,  the  rear  separation  will  enhance  the  pressure 
recovery  of  the  upstream  wake  and  the  enhanced  cavity  plateau 
pressure  will  in  turn  affect  the  downstream  separation  as  will 
be  explained  in  detail  later.  Whether  the  limiting  value  of 
reach  a unique  value  or  not,  when  either  the  step 
height  approaches  zero  or  the  cavity  length  approaches  infinity, 
can  only  be  concluded  by  further  tests. 
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Figure  44.  Effect  of  depth  to  length  ratio 
on  the  plateau  to  free-stream 
pressure  ratio. 
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4.1.6  Discussion  of  the  Separation  Pressure 

Through  a fairly  wide  survey  of  two-dimensional  data  on 
the  pressure  rise  to  separation  in  supersonic  and  high  tran- 
sonic flow,  it  has  been  reported  that  the  separation  pressure  is 
nearly  independent  of  Reynolds  number  at  sufficiently  large 
Reynolds  number  (Refs.  11,  12  and  13).  In  incompressible  flow, 
the  separation  pressure  coefficient  was  reported  to  be  a function 
of  Reynolds  number  if  this  is  sufficiently  low  (Refs.  14  and  15). 

In  the  present  experiments,  it  has  been  found  that  the  non- 
dimensional  separation  pressure  (Pg/p^)  has  a rather  strong 
dependence  on  Reynolds  number  (Figure  45) . It  should  be  cautioned 
however  that  there  are  small  variations  of  Mach  number  which 
accompany  the  changes  in  Reynolds  number  and  thus  make  it  difficult 
to  separate  clearly  the  exact  effect  of  Reynolds  number  on  sepa- 
ration at  lower  transonic  and  high  subsonic  speeds.  However, 

Figure  45  shows  clearly  a trend  of  the  non-dimensional  separation 
pressure  increasing  with  Reynolds  number  up  to  Re.  * 30  x 10 6 

and  this  is  followed  by  a region  of  near  independency.  This 
tendency  is  true  also  for  cavities  of  different  depth  to  length 
ratios . 

The  effect  of  the  free- stream  Mach  number  (Mot)  upon  the 
non-dimensional  separation  pressure  ratio  is  shown  in  Figure 
46.  It  is  found  that  Pg/p*,  also  increases  with  Mm.  The  dis- 
persion of  data  points  in  this  figure  is  probably  due  to  the 
effect  of  Reynolds  number.  Figure  46  also  shows  that  the  magni- 
tude of  the  non-dimensional  separation  pressure  will  be  higher 
if  the  flow  experiences  an  upstream  disturbance.  This  result  is 
equivalent  to  that  of  the  vortex  generator  performance  of  the 
rearward-facing  step  case  discussed  in  the  last  subsection.  The 
increase  of  the  non-dimensional  separation  pressure  ratio  with 
free-stream  Mach  number  for  the  flow  approaching  a forward-facing 
step  with  a natural  oncoming  boundary  layer  will  follow  the  same 
trend  as  in  supersonic  flow.  This  result  has  been  shown  in  Figure 
30,  as  discussed  previously. 
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Figure  46.  Effect  of  free-stream  Mach  number  on  the 
separation  to  free-stream  pressure  ratio. 
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Since  tine  upstream  disturbance  may  be  initiated  by  numerous 
sources,  it  is  sometimes  more  convenient  to  use  the  flow  para- 
meters at  certain  locations  downstream  of  the  rearward-facing 
step  as  a reference  value  to  predict  the  phenomenon  of  the  down- 
stream separation.  If  the  flow  properties  at  the  starting  point 
of  the  cavity  pressure  plateau  are  used  as  the  reference  upstream 
conditions,  then  the  smaller  the  magnitude  of  the  separation- 
plateau  pressure  ratio  (p  /p  ) the  more  difficult  it  becomes  for 

S p 

the  flow  downstream  to  separate  under  the  same  free-stream  con- 
ditions. This  result  is  shown  in  Figure  47  where  the  free-stream 
pressure  is  considered  as  the  plateau  pressure  for  the  forward- 
facing step  model. 

The  separation-plateau  pressure  ratio  will  increase  mono- 

tonically  with  free-stream  Mach  number  (Figure  47) , but  only  up 

to  Re.  = 30  x 10 ^ (Figure  48).  This  means  that  a continuous 
R 

increase  of  free-stream  Mach  number  can  continuously  enhance  the 
ability  of  the  flow  to  overcome  a downstream  separation.  But, 
with  respect  to  the  Reynolds  number,  there  exists  an  optimized 
value  which  can  suppress  the  flow  separation  with  the  highest 
efficiency.  The  effect  of  different  Reynolds  numbers  upon  flow 
separation  will  be  enhanced  as  the  free-stream  Mach  number  exceeds 
0.8  as  can  be  realized  from  examining  Figures  47  and  48.  Based 
on  the  results  of  this  study,  it  can  be  stated  that  this  con- 
clusion holds  at  least  for  the  range  0.6  < M <0.9  and  16  x 106 

< Re,  < 37  x 10 6 . 

LR 

The  separation  pressure  to  base  pressure  ratio  (p  /pfa) 

characterizes  the  pressure  energy  supplied  from  the  external  flow 

by  the  turbulent  mixing  process  to  the  viscous  layer  during  the 

journey  of  flow  along  the  cavity.  Figure  50  shows  that  increase 

of  Reynolds  number  can  increase  the  rate  of  energy  supply  in  the 

range  of  15  x 10 6 .<  Re,  < 30  x 106,  In  the  range  of 

30  x 10  £ Re^  x<  35  x 10  , the  rate  of  energy  supply  is  independent 
R 

of  the  Reynolds  number.  That  an  increase  of  free-stream  Mach 


95 


A EDC-TR  -77-1  03 


□ Forward-facing  step 
h = 1 in. , Lf=  74  in. 
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Figure  47 


Effect  of  free-stream  Mach  number  on  the 
separation  to  plateau  pressure  ratio. 


M =0.80 


Figure  50.  Effect  of  Reynolds  number  on  the  separation 
base  pressure  ratio. 
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number  can  increase  the  rate  of  energy  supply  is  shown  in  Figure 
51.  From  this  figure,  it  can  be  found  again  that  the  Reynolds 
number  effect  is  enhanced  as  the  free-stream  Mach  number  exceeds 
0.8. 

The  depth  to  length  ratio  of  the  cavity  is  another  parameter 
which  will  affect  the  wake-separation  process  for  flow  over  a 
cavity- like  model.  From  Figures  52  and  53,  it  can  be  realized 
that,  as  far  as  the  rear  separation  is  concerned,  there  exists 
an  optimized  depth  to  length  ratio  for  the  flow  over  a shallow- 
cavity  model  at  a given  Mach  number  and  Reynolds  number.  Usually, 
for  flow  over  a discontinuity  on  the  boundary,  there  will  be  some 
energy  loss  which  may  reduce  the  peak  pressure  ahead  of  the  down- 
stream step  to  a degree  lower  than  that  of  the  flow  approaching 
a forward- facing  step  without  the  upstream  disturbance.  But  an 
upstream  disturbance,  if  located  at  some  proper  distance  ahead 
of  the  rear  separation,  can  act  as  a vortex  generator  to  energize 
the  viscous  layer.  With  this  in  mind,  it  can  be  understood  how 
the  rear  separation  can  occur  at  a higher  pressure  condition 
(Figure  53).  For  the  flow  at  a Mach  number  of  0.6  and  a Reynolds 
number  at  15  x 106,  the  optimized  depth  to  length  ratio  (h/i)  is 
about  0.025  (corresponding  to  4 = 23.5  in,).  It  has  been  pointed 
out  earlier  that  the  flow  will  become  unsteady  and  the  cavity  is 
ready  to  open  if  the  depth  to  length  ratio  is  less  than  0.0634 
(£  = 15.75  in.).  Therefore,  in  this  case  the  optimized  depth  to 
length  ratio  of  the  cavity  is  one-third  larger  than  the  critical 
depth  to  length  ratio  under  the  same  free-stream  conditions. 

That  a disturbed  flow  can  sometimes  have  a peak  pressure  higher 
than  the  undisturbed  flow  can  be  explained  as  follows:  as  it 
turns  over  the  rearward-facing  step,  the  flow  speed  increases 
drastically  and  this  is  accompanied  by  a large  pressure  drop  and 
the  flow  becomes  highly  turbulent.  Due  to  this  high  turbulent 
mixing  phenomenon,  the  energy  of  the  external  flow  is  pumped 
into  the  boundary  layer  at  a very  high  rate,  The  surface 
pressure  in  the  reattachment  region  then  recovers  from  the  base 
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Figure  510  Effect  of  free-stream  Mach  number  on  the  separation 
to  base  pressure  ratio. 
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Figure  53. 


Effect  of  depth  to  length  ratio  on  the 
separation  to  free-stream  pressure  ratio. 
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pressure  to  the  plateau  pressure  with  a very  steep  gradient. 

The  plateau  pressure  always  attains  a level  higher  than  that  of 
the  free-stream  pressure.  In  addition  to  this  high  pressure 
level,  there  is  a continuous  inflow  of  high  energy  fluid  from 
the  external  region  at  a rate  higher  than  that  of  an  equilibrium 
flow.  This  results  in  the  total  energy  of  the  boundary  layer  in 
the  region  slightly  downstream  from  the  start  of  the  cavity  pres- 
sure plateau  being  higher  than  that  of  the  corresponding  flow 
approaching  the  downstream  step  in  the  absence  of  the  upstream 
disturbance.  Therefore,  the  disturbed  flow  can  overcome  the 
downstream  separation  easier  than  the  undisturbed  flow,  if  the 
forward- facing  step  is  located  immediately  downstream  of  this 
region.  Downstream  from  this  region,  the  rate  of  energy  in- 
flow slows  down  due  to  the  decrease  of  the  turbulence  level. 

The  shear  layer  velocity  profile  recovers  continuously  by 
absorbing  the  energy  from  the  inflow  of  fluid,  but  the  loss  due 
to  viscous  dissipation  remains.  Therefore,  the  net  inflow  of 
energy  which  can  benefit  the  surface  pressure  becomes  negative. 

As  the  flow  moves  further  downstream  to  a certain  distance,  the 
magnitude  of  the  surface  pressure  may  be  lower  than  that  for  a 
flow  without  the  upstream  disturbance.  Therefore,  for  a cavity 
of  that  length  or  longer,  the  peak  pressure  ahead  of  the  down- 
stream step  will  be  lower  than  that  of  the  undisturbed  flow 
approaching  a forward- facing  step.  On  the  other  hand,  if  the 
cavity  is  too  short,  the  quantity  of  the  incoming  energy  from  the 
external  region  will  not  be  enough  to  raise  the  energy  in  the 
viscous  layer  to  a level  higher  than  that  of  the  flow  without  an 
upstream  disturbance  during  the  limited  journey  of  the  flow  in 
the  cavity,  So,  the  peak  pressure  ahead  of  the  downstream  step 
will  also  be  smaller  than  that  for  a flow  without  the  upstream 
disturbance.  Only  when  the  depth  to  length  ratio  is  within  a 
proper  range  does  the  boundary-layer  flow  have  more  ability  to 
overcome  the  downstream  separation  than  does  the  flow  approaching 


104 


A E DC -TR  -77-1  03 


a forward- facing  step  with  an  undisturbed  natural  boundary  layer 
upstream. 

4.2  CORRELATION  OF  THE  SEPARATION  PRESSURE 

A purely  empirical  equation  for  the  prediction  of  the 
separation  pressure  ahead  of  a forward-facing  step  with  an  up- 
stream disturbance  has  been  established  as  follows: 


ReT  x 10 
lr 


~ 6 


1/4 


= 2 


(4.3) 


This  equation  is  plotted  in  Figure  54.  It  shows  good  agreement 

with  the  experimental  results.  The  maximum  error  is  about  2.5 

percent  and  occurred  at  the  depth  to  length  ratio  h/£  = 0.0425 

(£  = 23.5  in.).  This  is  the  ratio  which  results  in  the  maximum 

separation  pressure  for  a given  flow  condition  (see  Figure  53) . 

By  fixing  attention  on  the  data  points  of  flow  at  M =0.6  and 
6 00 
ReL  = 15  x 10  , it  can  be  found  that  the  trace  of  deviation  of 
R 

the  data  points  from  the  empirical  equation  has  the  same  trend  as 
that  shown  in  Figures  52  and  53.  Therefore,  it  can  be  reasonably 
concluded  that  the  error  of  this  equation  is  mostly  due  to  the 
neglect  of  the  depth  to  length  ratio  which  is  another  parameter 
that  can  affect  the  separation  pressure  as  pointed  out  in  the 
last  subsection. 

4.3  TURBULENT  MEAN-VELOCITY  MEASUREMENTS 
4.3.1  Mean  Velocity  Profiles 

The  total  pressure  distribution  across  the  shear  layer  and 
the  corresponding  static  pressure  were  used  to  yield  the  Mach 
number  and  the  velocity  distribution  at  different  stations  along 
models.  The  following  equations  were  used, 
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and. 
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1-1/2 


-1 


M 
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oo 


1 + 0.2  r M^1/2 
1 + 0,2  r M2  I 


(4.4) 


(4.5) 


where  r is  the  recovery  factor  which  was  assumed  to  be  related 
to  the  Prandtl  number  by  the  relation 


r 


(Pr) 


1/3 


In  this  experiment  emphasis  was  placed  on  the  study  of  the 
wake-separation  interaction  phenomenon  instead  of  the  vortex- 
mechanics,  as  was  done  in  most  of  the  published  literature.  After 
extensive  testing  on  different  depth  to  length  ratio  cavities, 
it  was  found  that  cavities  which  have  a depth  to  length  ratio  of 
0.0425  to  0.0250  (corresponding  to  Z - 23.5  and  50  inches  for 
h * 1 inch)  are  the  most  suitable  geometries  for  the  present 
study.  These  cavities  revealed  the  most  information  on  the  wake- 
separation  interaction  phenomenon.  They  are  long  enough  to  avoid 
the  dividing  streamline  jumping  across  the  entire  upstream  wake 
and  the  downstream  separation.  Therefore,  in  this  second  phase 
of  the  study,  cavities  of  40- , 27.5-  and  23, 5- inches  long  with  a 
1-inch  depth  were- used  as  the  principal  models.  The  stations  for 
the  velocity  profile  measurement  were  closely  spaced  along  these 
models.  Two  shorter  cavities  (h f%  = 0.0909  and  0.0634,  corres- 
pondong  to  Z = 11  and  15.75  inches)  were  also  used  to  study  the 
feature  of  cavity  flow  under  the  situation  when  the  dividing 
streaming  was  about  to  jump  across  the  cavity.  Only  at  one  or  two 
stations  were  velocity  profiles  measured  for  these  short  cavities. 
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Mach  number  versus  a vertical  distance  y is  plotted  in 
Figures  55  through  61,  The  continuous  variation  of  Mach  number 
along  the  flat  plate,  the  forward-facing  step,  the  40-,  27,5- 
and  23.5-inch  long  cavity  models  is  shown  in  Figures  62  through 
66,  respectively.  At  the  reattachment  points,  velocity  profile 
was  about  three  inches  upstream  of  the  separation  point.  At 
this  station  (3M/9y>w  f 0 was  expected,  A special  feature  of  the 
velocity  profile,  near  the  step-induced  separation  point,  needs 
attention.  The  profile  is  slightly  different  from  that  of  a 
natural  separation  by  a uniform  adverse  pressure  gradient  (in  the 
flow  direction).  The  velocity  profile,  close  to  a step-induced 
separation  point,  shows  an  abnormally  thin  portion  (velocity 
defect)  roughly  at  one  step  height  from  the  wall  surface.  This 
is  due  to  the  non-uniformity  in  traverse  static  pressure  variation. 
For  incompressible  and  supersonic  flow,  the  static  pressure 
distribution  ahead  of  the  step  increases  upward  from  the  wall 
and  reaches  a maximum  at  roughly  y/h  - 1 (Refs.  16  and  17). 

In  transonic  and  subsonic  flow,  this  phenomenon  is  also  expected. 
Alo^S  the  middle  portion  of  the  cavity,  the  velocity  profile 
continuously  becomes  fuller  downstream  from  the  reattachment 
point  and  starts  to  thin  again  (u/uro  decreases  for  a given  y/<5) 
when  it  is  near  the  downstream  separation  point.  The  recovery 
and  the  decay  of  these  velocity  profile  changes  are  rather 
gradual  as  compared  with  the  pressure  variation  along  the  model 
surface . 

The  dimensionless  velocity  profiles  (u/uj , at  the  middle 
station  between  the  reattachment  point  and  the  separation  point, 
for  flow  over  different  cavity  models  with  depth  to  length  ratios 
h/ £,  = 0.0909,  0.0634,  0.0425  and  0,0360  (i.e.,  i = 11,  15.75, 

23.5  and  27.75  inches)  are  plotted  against  y/<$  in  Figure  67. 

The  velocity  profile  for  flow  over  a flat  plate  at  the  station 
61  (indicating  the  distance  measured  from  the  leading  edge  in 
inches)  is  also  shown  in  Figure  67  for  comparison.  Data  for  a 
cavity  of  h/i,  = 0.0909  (£  = 11  in.)  have  shown  some  dispersion 
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Figure  58.  Velocity  profiles  on  the  27.5-inch 
shallow-cavity  model. 


12 


AEDC-TR-77-103 


M 

Figure  61.  Velocity  profile  on  the  11-inch 
shallow-cavity  model. 
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= 0.6,  ReL^  = 32.5  x 106 

Figure  62.  Flow  variation  along  the  flat-plate  model  (The 
Reynolds  number  is  calculated  for  the  distance 
to  the  forward-facing  step  in  Figure  63;  from 
AEDC-TR- 75-95) . 
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Me  = 0.6,  ReL^  = 32.5  x 106 

Figure  63.  Flow  variation  along  the  forward-facing  step 
model  (from  AEDC-TR-75-95) . 
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Figure  64.  Flow  variation  along  the  40-inch  cav 


ity  model 
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Figure  66.  Flow  variation  along  the  23. 5- inch  cavity  model 
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at  y/6  < 1.5,  which  is  due  to  the  unsteadiness  of  the  flow  as 
the  cavity  flow  becomes  nearly  open,  i.e.,  the  dividing  stream- 
line is  about  to  leave  the  floor. 

The  velocity  profile  in  Figure  67  gets  fuller  with  decreasing 

depth  to  length  ratio  until  h/*,  = 0.0425  (£  = 23.5  in,)  and  then 

becomes  thinner  again.  As  mentioned  in  the  last  section,  the 

rearward-facing  step  of  a cavity  of  h/ l **  0.0425  (S,  = 23.5  in.) 

has  the  best  efficiency  in  energizing  the  shear  layer  of  the  flow 

at  Mco  ~ 0-6  and  Re^  = 15  x 106.  Therefore,  the  flow  over  this 

R 

cavity  will  have  a higher  energy  and  a fuller  velocity  profile 
than  the  same  flow  over  the  27.5-inch  cavity  at  the  station 
14.5  inches  downstream  from  the  rearward- facing  step. 

Figure  68  shows  the  effect  of  an  upstream  disturbance  in- 
duced by  the  rearward-facing  step  of  a cavity  with  a depth  to 
length  ratio  h/ Z = 0.0425  C&  = 23.5  in.).  The  upstream  disturb- 
ance makes  the  velocity  profile  in  the  region  of  y/<5  i 0.4  fuller 
than  it  would  be  without  the  upstream  disturbance.  For  a depth 
to  length  ratio  of  0.036  (Z  = 27.5  in.),  the  profile  almost  co- 
incides with  the  upstream  disturbance  free  case--corresponding 
to  a naturally  developed  boundary- layer  flow  approaching  a 

forward- facing  step.  For  = 0.84  and  Re,  = 24.8  x 106,  the 

LR 

rate  of  decrease  of  the  kinetic  energy,  in  the  region  close  to  the 
floor  of  the  cavity  of  h/£  = 0.0425  (S,  =23.5  in.),  is  higher  than 
that  for  the  flow  at  M = 0.56  and  Re,  = 15.8  x 106  (Figure  69) 

Lr 

The  results  shown  in  Figures  67,  68,  and  69  imply  that  the 
recovery  of  boundary  layer  properties  after  any  upstream  dis- 
turbance, such  as  change  of  free- stream  Mach  number  and  approach 
Reynolds  number,  or  a discontinuity  in  the  boundary,  etc.,  all 
start  in  the  shear- layer  adjacent  to  the  wall  surface  and  spread 
upward  as  the  flow  moves  downstream, 
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4.3.2  Law  of  the  Wall  and  Law  of  the  Wake 


Very  near  the  wall  surface  the  pitot-probe  measurements 
become  difficult  and  the  results  are  doubtful  due  to  the  com- 
bined effects  of  wall  interference  and  local  flow  turbulence. 
Therefore,  in  the  present  study,  the  wall  and  wake  law  in  the 
form 


luy] 

- r 1 

T 

+ B + -u  K 

\ " J 

K [5 

C4.6) 


was  fitted  by  using  the  outer  region  data,  with  k and  B taken 
as  0,41  and  5,0  respectively.  The  law  of  the  wake  was 
represented  by, 


1 


cos 


(4.7) 


The  compressibility  effect  was  taken  into  account  by  use  of  Van 
Driest rs  generalized  velocity  u*  as  a substitute  for  u (Ref.  18), 
i.e,  , 


where 
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and  r is  the  recovery  factor,  u and  M denote  local  condition 

e e 

at  the  external  edge  of  the  boundary  layer, 

Making  these  substitutions,  Equation  4,6  becomes 
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The  displacement  thickness  can  be  expressed  as 


5* 

5 


C4.ll) 


By  combining  this  with  Equation  (4,6)  yields  the  result, 
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(4.12) 


To  use  this  equation  the  displacement  thickness  6*  was  found  by 
a graphical  method  from  a plot  of  the  corresponding  velocity 
profile  and  and  ? were  determined  through  an  iteration  pro- 
cess. The  slope  of  the  wall  and  wake  law  at  the  outer  edge  of 
the  boundary  layer  is , 
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V I y + * i 

L 


(4.13) 


and  is  different  from  the  actual  flow  condition.  In  order  to 
fit  the  curve  with  the  experimental  data,  the  value  of  <5  was 
determined  by  a method  of  trial  and  error.  The  fitting  region 
is  usually  limited  to  ( y / 6)  < 0.9  for  flow  with  a strong  wake 
or  even  (y/S)  < 0.75  with  a weak  wake. 

A typical  plot  of  the  wall  and  wake  law  for  flow  over  a 
flat  plate  is  shown  in  Figure  70.  The  data  points  near  the  wall 
show  a value  which  is  too  high.  This  is  probably  due  to  the 
deflection  of  the  streamlines  inducted  by  the  trailing  vortex  ' 
system  which  is  generated  by  the  circular  tip  of  the  pitot-probe. 
The  calculated  wake-strength  parameter  tT  is  0.72  which  is 
slightly  greater  than  the  value,  ¥ = 0.55,  predicted  by  Coles 
(Ref.  19)  based  on  Wieghardt  and  Tillmann’s  incompressible  flow 
data  (Ref.  20).  The  higher  magnitude  of  it  of  the  present  study 
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is  probably  due  to  the  effects  of  the  higher  surface  roughness 
of  the  model  used  here  as  compared  with  the  waxed-plywood  plate 
used  in  Wieghardt  and  Tillman*s  test  (Ref,  20),  The  slight 
discontinuity  of  the  joints  between  the  floor  blocks  and  between 
the  filling  blocks  may  also  result  in  an  increase  of 7,  However, 
the  tests  of  all  the  different  models  in  the  present  study  are 
carried  out  under  the  same  roughness  condition.  Therefore,  the 
results  of  this  flat-plate  model  are  still  adequate  for  use  as  a 
standard  for  the  discussion  of  the  influence  of  wake  and/or 
separation  on  the  velocity  distribution  in  both  the  transverse 
and  streamwise  directions. 

The  wall  and  wake  law  plots  for  flow  over  a forward-facing 
step  model  are  shown  in  Figure  71.  By  comparing  the  results  for 
the  flow  over  a flat  plate  and  a forward-facing  step  model 
(Figure  72),  it  is  found  that  at  station  61  the  two  profiles 
nearly  coincide.  Only  around  the  portion  where  3000  < (uxy/v) 

< 5000  does  the  latter  show  a very  small  hump.  This  hump  will 
become  more  obvious  as  the  profile  is  taken  closer  to  the  step 
(Figure  71).  It  is  certainly  due  to  the  upstream  effect  of  the 
downstream  separation.  Another  comment  will  be  mentioned  here, 
namely,  that  the  Mach  number  at  the  outer  edge  of  the  boundary 
layer  for  the  forward- facing  step  is  slightly  smaller  than  that 
of  the  flat  plate.  However,  the  free-stream  Mach  number  of  the 
flow  over  the  forward- facing  step  model  should  be  slightly  higher 
than  that  of  the  corresponding  downstream  Mach  number  at  the 
external  edge  of  the  boundary  layer.  Therefore,  it  is  believed 
that  the  free-stream  Mach  number  for  these  two  profiles  is  nearly 
equal . 

For  flow  over  the  40-  and  27.5-inch  cavities,  the  plots  of 
the  wall  and  wake  law  are  shown  in  Figures  73  and  74,  respectively. 
After  examining  these  figures,  several  points  can  be  made  as 
follows : 
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Figure  71.  Wall  and  wake  law  plots  for  flow  over  the 
forward- facing  step  model. 
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Figure  73 „ Wall  and  wake  law  plots  for  flow  over  the 
40-inch  cavity  model. 
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Figure  74.  Wall  and  wake  law  plots  for  flow  over 
the  27.5-inch  cavity  model. 
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(a)  The  velocity  profiles  recover  continuously  from  the 
effects  of  the  rearward-facing  step  after  passing  the  reattach- 
ment point  in  the  40-inch  cavity  case.  The  upstream  effect  of 
the  downstream  separation  will  overwhelm  the  relaxing  process 
after  the  flow  has  reached  the  station  70  (Figure  73) , For  flow 
over  the  27.5-inch  cavity,  the  velocity  profiles  at  fewer 
stations  were  measured;  hence  no  such  conclusion  can  be  made. 

(b)  At  the  location  21  inches  downstream  from  the  reattach- 
ment point  (station  61),  in  the  40-inch  cavity,  the  outermost 
portion  of  the  wall-  and  wake-law  profile  coincides  closely  with 
that  of  the  upstream-disturbance-free  case,  as  is  clearly  illus- 
trated in  Figure  75.  In  the  region  0.07  < (y/5)  <0.3  the 
forward- facing  step  profile  has  a substantial  hump.  Since  the 
external  flow  velocity  of  the  two  flows  is  equal  and  the  frictional 
felocity  (u^)  of  the  cavity  flow  is  18.8  ft/sec  while  it  is 

20.6  ft/sec  for  the  forward-facing  step,  it  is  concluded  that  this 
hump  denotes  the  kinetic  energy  of  the  cavity  flow  near  the  floor 
must  be  smaller.  Furthermore,  the  surface  pressure  is  14.6  psia 
for  the  cavity  flow  and  is  15.2  psia  for  the  forward-facing  step. 
Therefore,  the  total  energy  of  the  cavity  is  smaller  and  this 
results  in  a lower  peak  pressure  ahead  of  the  downstream  step  as 
was  pointed  out  in  Section  4.1,0. 

(c)  At  stations  near  either  the  reattachment  or  the  sepa- 
ration point,  the  agreement  with  the  law  of  the  wall  is  fairly 
poor.  Vasanta,  Ram  and  Wauschkuhn  (Ref.  21)  assumed  that  at 
different  stations  downstream  from  the  reattachment  point  the  law 
of  the  wall  should  be  fitted  by  different  slopes.  Schofield  (Ref. 
22)  thought  that  this  change  of  slope  is  attributed  to  this  "non- 
equilibrium of  the  wall  turbulence,"  Coles  (Ref.  23)  said  that 
near  the  wall  data  could  be  wrong  but  the  formula  is  correct. 

The  present  experiments  show  a tendency  for  the  slope  of  the  law 
of  the  wall  at  the  reattachment  point  to  be  smaller  than  that  of 
the  flat  plate.  The  slope  increases  as  the  flow  moves  downstream 
and  finally  becomes  greater  than  the  flat  plate  case  as  the  flow 
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approaches  the  separation  point  (Figures  73  and  74) , The  data 
points  in  the  wall  region  more  nearly  fall  along  a straight  line 
than  along  a sinuous  shape  as  was  asserted  by  Schofield  (Ref.  22) 
based  on  Bradshaw  and  Wong's  experiments  (Ref.  24).  However,  the 
flow  measurement  using  an  impact  probe  near  the  wall  with  a high 
turbulence  usually  resulted  in  a higher  degree  of  wall  inter- 
ference and  thus  yielded  data  of  relatively  low  reliability.  Un- 
fortunately, all  existing  experiments  rely  on  such  a measuring 
technique;  hence  all  the  different  explanations  concerning  the 
agreements  with  the  law  of  the  wall  are  plausible,  but  open  to 
question.  Before  the  data  can  be  obtained  by  an  interference- 
free  type  of  measurement,  no  definite  conclusion  can  be  suggested 
here . 

(d)  The  law  of  the  wake  profile  shows  a hump  at  station 
61  in  the  40- inch  cavity  as  mentioned  in  (b)  above.  It  dis- 
appears at  the  station  66,  and  shows  up  again,  more  strongly,  at 
station  70.  In  the  27.5-inch  cavity  flow,  there  is  also  a 
significant  hump  on  the  wake  profile  at  the  station  70.  Since 
it  has  been  found  that  the  static  pressure  ahead  of  the  down- 
stream step  increases  transversely  from  the  floor  and  reaches 
a maximum  at  one  step  height  as  in  a supersonic  or  an  incom- 
pressible flow  passing  a shallow  cavity,  the  transverse  pressure 
gradient  becomes  larger  at  the  stations  closer  to  the  forward- 
facing step  (Ref.  16).  Although  there  has  been  no  measured 
information  about  this  for  subsonic  and  transonic  flow,  it  is 
expected  that  the  tendency  must  be  the  same,  This  transverse 
static  pressure  -gradient  will  produce  a change  in  the  velocity 
profile,  The  result  is  that  the  velocity  profiles  near  the 
downstream  step  show  an  unusual  velocity  defect  which  occurs 
between  a height  of  one  step  height  to  the  layer  near  the  floor. 
This  unusual  defect  region  forms  the  hump  on  the  wake-law  plot. 
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4.3.3  The  Wake  Strength  Parameter  it  and  the  Velocity  Profile 
Defect  Parameter  G 

The  law  of  the  wake  is  characterized  by  a universal  function 
a)  (y/6)  for  all  two-dimensional  turbulent  boundary  layers  and  a 
wake  strength  parameter  tT  which  specifies  different  strengths  of 
the  wake  component.  Coles  found  that  n is  very  nearly  0.55  for 
flow  at  constant  pressure  (Ref.  19)  but  is  different  under  dif- 
ferent non-equilibrium  flow  conditions  (Ref.  23).  tt  increases 
for  the  flow  over  an  adverse  pressure  gradient  field  and  decreases 
with  a favorable  pressure  gradient.  Its  magnitude  can  be  either 
positive  or  negative.  At  the  separation  or  reattachment  point, 
rF  approaches  infinity.  The  results  for  tt  from  the  present 
experiments  are  plotted  in  Figure  76. 

For  flow  over  a flat  plate,  Prandtl  derived  an  equation  for 
the  estimation  of  the  skin  friction  (Ref.  25),  i.e., 

Cf  = 0.074  (Rex)~1/5;  5 x 105  < Re^  < 107,  Re6  < 105 (4 . 14) 

By  substitution  and  rearrangement  from  the  above  equation,  and 
by  assuming  that  (5*/ 6)  = (1/8)  and  that  k = 0,41,  the  expres- 
sion for  tt  becomes , 


— uco  n 1 

{T  - 1 - °-3  CRex)U-L  -1  (4.15) 

This  equation  shows  that  it  for  the  flow  over  a flat  plate 

is  a function  of  Reynolds  number  only,  Intuitively,  tt  should 

depend  also  on  the  surface  roughness  which  will  show  up  as  a 

component  of  skin  friction.  By  observation  of  the  above  two 

equations,  it  is  understood  that  the  roughness  effect  is  embedded 

in  the  Reynolds  number,  The  Reynolds  numbers  based  on  length 

(Re  ) and  the  boundary  layer  thickness  (Re J for  the  present 
X 6 _ 

experiments  are  about  32.5  x 10°  and  6 x 10  , respectively,  for 
flow  at  Mro  = 0.6.  These  values  are  slightly  higher  than  the 
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magnitudes  specified  in  the  above  equation.  However,  the 
magnitude  of  the  wake  strength  parameter,  It = 0,69,  determined 
by  the  equation,  is  still  very  close  to  the  measured  data  of 
0.72.  This  result  is  much  better  than  that  obtained  from  the 
Schultz-Grunow  skin  friction  formula  fRef,  26) , i.e., 

Cf  - 0.427  Clog  Rex  - 0,407)~2,64  C4.16) 

i =0,11  (log  Re^  - 0.407) 1,32  - 1 

which  gives  a rather  low  value  of  ¥ CO. 46). 

Clauser  (Ref.  27)  introduced  a universal  thickness  para- 
meter, 


4 


u - u 

00 

U 


dy 


(4.17) 


for  turbulent  boundary- layer  data  analysis,  which  will  make  the 
velocity  defect  plot,  Cu  - uoo)/uT  versus  (y/A)  , independent  of 
Reynolds  number  and  the  surface  roughness.  By  substitution  of  ue 
for  u^  in  the  above  equation  and  by  the  application  of  the  wall 
and  wake  law  with  the  boundary  condition  u = ue  as  y ■+  5,  yields, 

" = k | - 1 (4.18) 


This  equation  shows  that  ¥ is  linearly  related  to  (.4/6),  where  6 
contains  the  mixed  response  of  Reynolds  number,  surface  roughness 
and  any  other  existing  upstream  disturbances.  The  plot  of  ¥ 
versus  (A/ 6)  is  shown  in  Figure  77. 

In  addition  to  A,  Clauser  introduced  a velocity  profile 
defect  parameter  G to  characterize  the  various  members  of  the 
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universal  family  of  turbulent  boundary- layer  velocity  profiles 
in  the  following  fashion,  i.e,, 


G = / 

o 


u - u 


1 UT  i 


dy  / / 

0 


°°/u „ - U ' 


u 


dy 


= / 
o 


/ um  - u 

; 

00 

1 d f£! 

u 

J ‘lAI 

T ( 

(4,19) 


Similarly, 


G can  be  expressed  as 

A- 


, ll2 

1 + COS  T7  (^J  | ! d 


(4.20) 


After  integration  (Refs.  28  and  29),  the  last  equation  yields, 


G = | (11.89  + 24,85  tr  + 8.92  tt  2)  (4.21) 

The  plot  of  G versus  (A/6)  is  also  shown  in  Figure  77.  This 
equation  gives  G = 7.47  corresponding  to  tT  = 0.55  as  compared  to 
G = 6.1  obtained  by  Clauser  (Ref,  27). 

4.3,4  The  Velocity  Defect  Law 

In  the  preceding  subsection,  it  was  shown  that  the  law  of 
the  wall  and  the  law  of  the  wake  can  yield. 


r 

\Z\  . ZL 

i 1 

j 1 + cos 

IT 

6 K 1 

1 

6 1 

1 

L. 

■ 


(4.22) 


This  equation  is  known  as  the  velocity  defect  law,  The  plot  of 
the  defect  law  for  the  flow  over  the  40-inch  cavity  is  shown  in 
Figure  78.  It  reveals  the  same  trends  as  those  exhibited  by  the 
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Figure  78.  Velocity  defect  plot  for  flow  over  the 
40-inch  cavity  model. 
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wall  and  wake  law,  namely,  that  the  upstream  effect  of  the  down- 
stream separation  does  not  overwhelm  the  relaxation  process  until 
station  70,  Figure  79  is  a velocity  defect  plot  at  the  station 
70  for  the  flow  over  different  models--the  40-  and  27,5-inch 
cavities,  the  forward- facing  step  and  the  flat  plate,  It  shows 
the  phenomena  of  the  downstream  effect  on  the  upstream  disturb- 
ance and  the  upstream  effect  on  the  downstream  separation.  The 
inner  portion  of  the  velocity  defect  profile  for  the  flow  over 
a forward-facing  step  model  is  slightly  fuller  than  that  of  the 
flow  over  the  27,5-inch  cavity  model  but  thinner  than  that  for 
the  flow  over  a 40- inch  cavity  model.  Unlike  the  relaxation  of 
the  base  pressure  which  will  slightly  overshoot  almost  immediately 
downstream  of  the  reattachment  point  and  then  decreases  slightly 
as  the  flow  moves  further  downstream,  the  velocity  profile  will 
relax  at  a slower  rate,  At  a station  about  18  inches  downstream 
from  the  reattachment  point  (i,e.,  station  70  in  the  27.5-inch 
shallow  cavity) , the  inner  portion  of  the  velocity  profile  has 
recovered  to  the  state  of  the  upstream  disturbance-free  case 
(i.e.,  flow  over  a forward-facing  step  model).  While  at  the 
same  station,  but  with  a disturbance  at  about  30  inches  upstream 
(i . e . , flow  over  40-inch  shallow  cavity),  the  inner  portion  of 
the  profile  becomes  even  fuller  than  that  without  an  upstream 
disturbance  and  almost  coincides  with  the  profile  for  the  flow 
over  a flat  plate.  This  is  perhaps  due  to  the  turbulent  mixing 
process  which  can  transport  energy  from  the  external  flow  into 
the  inner  shear  region. 

It  would  be  of  interest  to  determine,  when  the  cavity  length 
is  further  increased,  whether  or  not  the  turbulent  mixing  pro- 
cess can  cause  the  velocity  profile  at  certain  downstream  stations 
to  become  fuller  than  the  flat  plate  case  followed  by  a subse- 
quent thinning.  As  far  as  the  outer  portion  of  the  viscous  layer 
is  concerned,  the  velocity  defect  profile  coincides  with  that 
of  the  flow  over  a flat  plate  if  there  is  no  upstream  disturbance, 
The  closer  the  upstream  disturbance  is,  the  thinner  will  be  the 
outer  portion  of  the  velocity  defect  profile. 
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4.4  BOUNDARY -LAYER  THICKNESS  PARAMETER 

The  boundary- layer  thickness  ($)  development  along  different 
models  (flat  plate,  forward-facing  step,  27.5-  and  40-inch 
cavities)  for  flow  at  a free-stream  Mach  number  = 0,6  and 
unit  Reynolds  number  (u^/v)  = 4.5  x 10^  are  shown  in  Figure 
80.  These  values  were  determined  by  the  rule  u « 0.99  ue  at 
y = <5 - The  boundary-layer  thickness  defined  by  the  law  of  the 
wall  and  the  law  of  the  wake,  known  as  the  Coles  boundary- layer 
thickness  ( 6 c ) , are  also  shown  in  this  figure.  When  the  wake 
component  is  small,  the  difference  between  5-  and  <$'c  is  significant. 

The  displacement  thickness  ($*)  was  determined  by  a 
graphical  method  from  the  plot  of  the  corresponding  velocity 
profile  according  to  the  definition 


6* 


CO 


The  displacement  thickness  based  on  Coles  boundary  layer  thick- 
ness ($c)  was  also  determined.  The  results  obtained  for  6* 
are  shown  in  Figure  81. 

The  momentum  thickness  is  calculated  as  follows . 


e 


= 7c~Ji  (l  - gT5^ T a 
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where  Cf  is  determined  by  the  law  of  the  wall  and  the  law  of  the 
wake.  The  momentum  thickness  determined  by  this  equation  should 
be  called  Coles  momentum  thickness  (8  ) , The  results  determined 
for  Coles  momentum  thickness  are  shown  in  Figure  82. 

4.5  WALL  SHEAR  STRESS  DISTRIBUTION 

The  skin  friction  coefficient  (Cf)  was  determined  from  the 
result  of  the  law  of  the  wall  and  the  law  of  the  wake,  The 
direct  measurement  of  the  skin  friction  coefficient  was  just 
beginning  at  this  reporting  time,  Therefore,  the  comparison 
was  done  through  the  wall  and  wake  law.  The  skin  friction 
variation  along  the  model  surface  is  shown  in  Figure  83.  In  the 
region  of  60  < x < 70  inches,  the  skin  friction  coefficient  for 
the  flow  over  the  40- inch  cavity  model  nearly  coincides  with  that 
for  the  flow  over  the  flat  plate  model  and  this  corresponds  to 
the  equilibrium  condition.  In  contrast,  the  plot  of  the  skin 
friction  coefficient  for  the  flow  over  the  27,5-inch  cavity 
always  shows  a lower  value  than  that  for  the  corresponding 
equilibrium  flow. 

The  plot  of  skin  friction  coefficient  versus  Reynolds 

JL 

number  based  on  Coles  displacement  thickness  (fi")is  shown 

c 

in  Figure  84.  The  empirical  relations  are  formulated  as 
follows ; 

For  the  flat  plate, 

C.  = 0.0189  (Re  *)  ~0-174 

c 

For  the  forward-facing  step, 

C,  = 2.38  (Re  *)'0’667 

r o 

c 
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Figure  82.  Development  of  Coles'  momentum  thickness 
along  different  models. 
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Figure  84.  Variation  of  skin  friction  coefficient 
with  a Reynolds  number  based  on  Coles' 
displacement  thickness. 
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Figure  85.  Variation  of  skin  friction  coefficient 
with  a Reynolds  number  based  on  Coles' 
momentum  thickness. 
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TABLE  III 

DIMENSIONS  OF  SKIN  FRICTION  BALANCE 
(See  Figure  86  for  Notation) 


Area  of  Circular  Floating  Element 
Width  of  Gap 
Length  L 
Length  Lj 


0.490  in2 
0.005  in 
9.25  in 
3,475  in 


Length  6 


0,625  in 
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For  the  40-inch  shallow  cavity, 

Cf  = 4,83  x 101Q  (Re^*)  ^*693,  41  < x < 51  inches 

y c 

Cf  = 1.13  CRe^*)  41  < x < 66  inches 

c 

Fox  the  27.5-inch  shallow  cavity, 

= 4,76  x 10^  (Re^*)  ^•^1)  41  < x < 61  inches 

c 

= 1.12  (Re^.*)  ^-567^  61  < x < 66  inches 

c 

Since  some  of  these  formulas  were  determined  by  only  two 
data  points,  the  reliability  is  somewhat  doubtful.  This  will 
be  checked  against  our  own  direct  force  measurements  soon. 

Figure  85  shows  the  relation  between  the  skin  friction  co- 
efficient  and  the  Reynolds  number  based  on  Coles  momentum  thick- 
ness, Reda's  (Ref-  30)  direct  measurement  of  the  surface  shear 
for  flow  over  a smooth  flat  plate  is  included  for  comparison. 

The  present  measurement  is  about  20  percent  higher  than  Redars 
result.  This  may  be  mainly  due  to  different  degrees  in  the 
surface  roughness  of  the  models  used. 

Skin  friction  coefficient  measurements  were  just  beginning 
at  the  end  of  this  reporting  period,  Therefore,  no  results  can 
be  presented.  However,  a complete  description  of  the  operation 
of  the  skin  friction  balance  will  be  given. 

The  operation  of  the  skin  friction  balance  is  explained  in 
detail  in  Figure  86.  Important  dimensions  are  given  in  Table  III, 
The  skin  friction  acting  on  the  area  of  the  floating  element 
causes  the  support  rod  to  rotate  about  point  A.  This  motion  is 
sensed  by  the  motion  of  the  LVDT  core,  at  point  E,  with  respect 
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to  the  LVDT  coil.  A pair  of  flexures.  Figure  87,  are  provided 
at  the  pivot  point.  The  flexures,  made  from  a single  piece  of 
stainless  steel,  consist  of  a center  axle  and  four  ribs  which 
radiate  from  the  axle  to  an  outer  rim.  The  rim  is  held  stationary 
in  the  balance  housing  while  the  axle  is  firmly  attached  to  the 
floating  element  arm,  When  the  arm  is  rotating  the  flexure  ribs 
deform  giving  rotational  stiffness  and  provide  a frictional 
pivot.  The  torque  provided  by  the  flexures  is  indicated  by  M 
in  Figure  86.  Damping  is  also  provided  in  the  flexure  to  lessen 
the  effects  of  vibrations  (see  Figure  87). 

A counterweight  C (Figure  86)  is  provided  to  place  the 
center  of  gravity  of  the  balance  arm  at  the  pivot  point  A.  This 
also  Bias  the  effect  of  lessening  the  sensitivity  of  the  balance 
measurement  to  tunnel  vibrations, 

Upon  the  application  of  a skin  friction  force,  F^,  the 
balance  arm  will  rotate  about  A,  displacing  the  LVDT  core  B. 

The  LVDT  signal  will  rotate  the  motor  causing  a linear  extension 
of  the  coil  spring.  The  elongated  spring  will  create  a force  F2 
which  will  increase  until  the  moment  of  the  skin  friction  force 
(F^  L)  will  equal  the  moment  of  the  spring  force  (F2b) . The 
total  rotation  of  the  motor  shaft,  as  measured  by  the  potentio- 
meter, is  then  proportional  to  the  force  F-^. 

The  electronics  which  control  the  sevomechanism  are  shown 
in  Figure  88.  The  LVDT  coil  is  electrically  excited  by  a 20  kHz 
signal  from  a Schaevitz  CAS- 200  signal  conditioner.  Movement  of 
the  LVDT  core  from  the  null  position  results  in  a coil  output 
signal  of  the  same  frequency  but  different  phase  from  the 
excitation  signal.  The  signal  conditioner  produces  a signed-DC 
output  proportional  to  the  AC  input  from  the  coil.  The  sign  of 
the  output  (plus  or  minus)  depends  upon  the  phase  of  the  coil 
output  with  respect  to  the  excitation  signal,  one  sign  for  each 
of  the  directions  of  motion  of  the  core  from  the  null  position. 

The  DC  signal  from  the  Schaevitz  signal  conditioner  is  fed 
to  the  DC  motor  controller,  The  motor  controller  consists  of 
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a high  current  operational  amplifier  (1  amp)  (figure  89) , The 
amplified  DC  signal  then  operates  the  DC  motor  (155  rpm,  60  in 
oz  maximum  torque) , The  speed  of  the  motor  is  controlled  hy  the 
magnitude  of  the  input  voltage,  which,  in  turn,  depends  upon  the 
distance  of  the  LVDT  core  from  null,  As  the  core  approaches  the 
null  position,  the  motor  will  slow  and  come  to  a stop,  if  enough 
damping  has  Been  provided  in  the  system  to  prevent  overshoot. 

The  skin  friction  balance  was  calibrated  by  hanging  small 
weights  from  a flexible  filament  which  was  attached  to  the 
floating  element  by  tape  and  suspended  over  a pivot  which  has 
jeweled  bearings  (Figure  90).  The  pivot  was  attached  to  the  end 
of  the  ’’dummy"  block  which  contains  the  floating  element;  the 
adjacent  floor  block  was  removed  to  facilitate  this  operation. 
Typical  calibration  curves  are  shown  in  Figure  91. 

A number  of  modifications  were  made  to  the  balance  before 
it  would  operate  in  the  UTSI  Transonic  Wind  Tunnel.  In  addition 
to  modifying  the  housing  to  attack  the  balance  to  the  UTSI 
tunnel,  modifications  also  were  made  to  the  mechanical  and 
electrical  parts  of  the  feed-back  loop. 

Two  aspects  of  the  operating  environment  in  the  UTSI  tunnel 
were  different  from  those  which  exist  in  the  NOL  tunnel,  namely, 
the  magnitude  of  the  skin  friction  force  to  be  measured  and  the 
vibrational  environment.  The  balance  was  originally  designed  to 
measure  skin  friction  forces  in  the  range  0.03  to  1 gin  while  the 
forces  expected  in  the  UTSI  tunnel  are  1 to  15  gm.  The  average 
UTSI  tunnel  acceleration  is  4 g's,  rms . Therefore,  a stiffer 
restoring  spring  and  additional  damping  were  required,  After 
replacing  a broken  flexure  and  the  damping  fluid  satisfactory 
performance  was  obtained  with  a very  stiff  restoring  spring, 

Two  motor  control  systems  were  examined.  The  first  motor 
controller  used,  manufactured  by  AST  Servo  Systems,  Inc,,  would 
not  operate  in  a stable  manner  and  was  replaced  by  a high  current 
operational  amplifier  (Figure  89),  Electronic  RC  damping 
circuits  were  not  required  to  stabilize  the  servo  circuit. 
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20  kft  pot 


Operational  Ampl. 


Figure  89,  Motor  controller. 
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Figure  91a.  Calibration  curve  of  ski 
meter  mounted  in  UTSI  Tr 
Wind  Tunnel. 
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The  signal  from  Che  potentiometer  which  was  geared  to  the 
servomotor  (Figure  10)  was  amplified  and  recorded  on  an 
oscillograph . 


5.0  THEORETICAL  CONSIDERATIONS 

In  order  to  correlate  the  experimentally  measured  pressure 
data  for  flows  past  cavities  at  different  Mach  numbers,  the 
following  theoretical  method  has  been  developed, 

y + 


> x 


In  a purely  subsonic  flow,  the  pressure  coefficient,  Cp  (x)  , 
for  flow  past  a cavity  as  shown  in  the  sketch  above,  based  on 
subsonic  small  disturbance  theory,  can  be  written  as 

£ I — ill 

-1  C-x 

We  have  developed  a method  (Ref,  1)  by  which  the  inverse  solu- 
tion for  the  above  equation  can  be  obtained,  that  is,  we  can 
obtain  an  expression  for  the  shape  of  an  inviscid  stream  line  on 
the  cavity,  the  flow  past  which  will  produce  an  experimentally 
measured  C^  - distribution  based  on  inviscid  subsonic  small  dis- 
turbance theory.  Such  an  expression  can  be  written  as 


Cp  (x) 
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where  r(x) 


{ (1  - x )/  ( 1 + x) 


This  expression  ought  to  give  a reasonably  good  approximation 
for  the  outer  inviscid  flow  stream  line  under  the  assumption 
that  the  normal  pressure  gradient  in  the  viscous  layer  is 
negligible.  Figure  92  represents  a typical  inviscid  stream 
line  shape  along  the  cavity  wall  which  corresponds  to  a 
measured  C (x). 


Based  on  a conjecture  that  the  thin  viscous  layer  thickness 
may  not  change  appreciably  when  we  increase  the  Mach  number  from 
a pure  subsonic  level  to  a transonic  level  (with  M /l,  still) 
we  have  used  the  outer  inviscid  stream  line  computed  as  indicated 
earlier  for  a pure  subsonic  flow,  to  compute  C (x)  at  a higher 
Mach  number  using  inviscid  transonic  small  disturbance  equation. 
The  small  disturbance  potential  v in  transonic  flow  is  governed 
by 


(K  " vx)  3xx  + = 0 * y = 5*y 

K = (1  - mVs1 

00 

We  have  modified  and  adapted  a computer  program  developed  by 
Murraan  and  Krupp  for  solving  such  an  equation  by  a mixed  finite 
difference  scheme  for  airfoil  flows,  to  solve  the  flow'  past  a 
cavity. 

An  experimental  C distribution  for  M = 0.58  and  cavitv 

P ca  J 

length  23.5"  has  been  utilized  to  compute  the  inviscid  stream 
line  along  the  cavity  by  the  inverse  method  for  pure  subsonic 
flows  outlined  earlier.  Using  this  inviscid  stream  line,  a 
transonic  inviscid  flow  computation  has  been  carried  out  for 
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Figure  92.  Computed  inviscid  streamline  corresponding  to 
the  experimentally  measured  C,,  distribution 

at  M = 0.58.  p 
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M = 0.84.  The  computed  C distribution  agrees  rather  well 

oo  p 

with  the  measured  C - distribution  at  M = 0.84  (Fig.  93). 

P 00 

This  shows  that  the  viscous  flow  is  limited  to  the  same  region 
at  both  Mach  numbers,  thus  the  inviscid  stream  line  outside 
the  boundary  layer  does  not  seem  to  change  much  for  Mach  numbers 
from  0.58  to  0.84. 


It  is  well  known  (Ref.  31)  that  the  inviscid  transonic 
equation  and  Rankine-Hugoniot  conditions  lead  to  inconsistencies 
in  the  region  where  the  shock  terminating  an  embedded  super- 
sonic pocket  touched  the  curved  surface  of  an  airfoil,  and  the 
computed  surface  pressure  jump  is  more  than  what  is  observed 
experimentally.  The  viscous-transonic  equation,  which  includes 
the  effects  of  compressive  viscosity  on  the  structure  and  thus 
location  of  the  shock,  seeks  to  overcome  these  shortcomings. 

It  is  a small  disturbance  approximation  to  the  Navier-Stokes 
equations  and  includes  the  effects  of  compressive  viscosity  on 
non- Rankine-Hugoniot  regions,  but  does  not  include  the  boundary 
layer  effects.  It  can  be  written  as 
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viscous  transonic  similarity  parameter 


This  equation  is  parabolic  everywhere,  unlike  the  inviscid 
transonic  equation  which  is  mixed  elliptic-hyperbolic,  but  it 
is  not  an  initial  value  problem  in  any  one  direction.  The 
boundary  conditions  are  prescribed  all  around  the  boundary  like 
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an  elliptic  problem.  We  have  adopted  a finite  difference 
scheme  which  is  similar  to  the  "box  scheme"  developed  by  Keller 
(Ref.  32).  The  non-linear  algebraic  equations  are  solved 
iteratively  by  a horizontal  line  relaxation  technique.  The 
results  are  very  promising.  The  details  of  the  method  of 
solution  are  outlined  in  the  Appendix. 

We  have  performed  test  calculations  for  a thin  parabolic 
symmetric  airfoil,  For  Ky  = 0.01  and  K.^  = 0.725,  starting  with 
zero  initial  conditions,  Figure  94  shows  the  development  of 
the  C distribution  on  the  airfoil  surface  at  100  sweep  intervals 
in  the  relaxation  method.  Figure  95  shows  the  converged  Cp 
after  1600  sweeps.  Several  runs  have  been  made  to  analyze  the 
influence  of  the  viscous  transonic  similarity  parameters  as  it 
is  reduced  to  smaller  and  smaller  values.  It  is  expected 
theoretically  that  the  results  ought  to  approach  those  of  the 
inviscid  flow  computations . It  has  been  found  that  when  the 
viscous  similarity  parameter  Kv  is  reduced  from  0.01  to  0.005 
and  to  0.001  for  a subcritical  case  with  Ki  = 1.67,  the  Cp  -distri- 
bution is  essentially  unchanged  and  it  agrees  with  known  com- 
putational and  experimental  results  (Figs.  96  and  97).  For  a 
super-critical  case  with  Y,i  = 1.0,  the  Cp  - distribution  has 
been  shown  to  undergo  little  change  as  Ky  is  reduced  from  0.01 
to  0.005  except  in  the  region  of  flow  re-expansion  behind  the 
shock  (Figs.  98  and  99).  The  method  has  also  been  extended  to 
compute  the  flow  with  a separation  wake  just  behind  the  shock 
terminating  the  sonic  rocket.  With  Ky  = 0.01  and  K±  = 1.00, 

Figures  100,  101  and  102  show  the  Cp  - distribution  with  or  with- 
out separation  wakes.  Figure  102  indicates  the  forward  movement 
of  the  shock  because  of  the  separation  and  the  significant  in- 
fluence on  the  flow  field  way  ahead  of  the  shock. 

The  feasibility  of  the  use  of  viscous  transonic  equation  for 
computing  invisid  transonic  flows  is  established.  The  method  of 
solution  has  some  advantages  over  the  method  of  mixed  difference 
schemes  for  inviscid  transonic  equation.  Further  applications  of 
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Figure  102.  Cp  - distribution  as  a parabolic  airfoil 
with  separation  wake  behind  the  shock  and 
the  sonic  pocket  - Kv  = 0.01,  K-j  = 1.0 
(14%  airfoil  at  M = 0.9). 
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the  method  and  more  refinements  of  the  method  itself  are  possible 
and  will  be  pursued. 


6.0  CONCLUSION 

The  second  phase  of  a program  to  study  the  wake/ separation 
interaction  flow  field  by  employing  several  shallow-cavity  models 
has  been  accomplished.  The  entire  program  is  intended  to  shed 
some  light  on  the  upstream  disturbance  effect  on  transonic  flow 
suffering  a significant  adverse  pressure  gradient  downstream. 

Both  experimental  and  theoretical  work  has  been  incorporated  in 
the  study.  It  is  hoped  that  through  this  program,  the  scaling 
effect  as  well  as  the  mechanism  by  which  the  flow  is  rehabilitated 
and  re-energized  can  be  better  understood. 

The  present  experiments  were  carried  out  with  a natural 
boundary  layer  developed  from  the  model  leading  edge.  The  extra- 
ordinarily long  model  results  in  a range  of  high  Reynolds  number 
which  can  cover  that  of  actual  flight  conditions.  The  upstream 
disturbance  was  simulated  by  a rearward-facing  step.  The  strength 
of  the  upstream  disturbance  was  controlled  by  varying  the  free- 
stream  Mach  number  and  Reynolds  number,  as  well  as  the  location 
of  the  steps. 

It  was  found  that  the  upstream  disturbance  can  either 
suppress  or  enhance  the  downstream  separation  depending  on  the 
magnitude  of  the  disturbance,  the  distance  between  them  and  the 
Reynolds  number  and  the  Mach  number  of  the  oncoming  flow.  The 
variations  in  the  base  pressure,  the  reattachment  pressure,  the 
plateau  pressure,  the  re-separation  pressure,  as  well  as  the 
peak  (i.e.,  maximum)  pressure  within  the  shallow-cavity  flow  by 
varying  various  parameters  have  been  obtained.  Among  these 
parameters,  it  was  found  that  the  flow  responded  more  acutely 
to  the  depth/ length  ratio  and  the  Reynolds  .number . 
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Based  on  the  present  experiments,  it  was  found  that  the 
Reynolds  number  influences  the  flow  pattern  strongly  initially, 
however,  it  then  leveled  off  at  a certain  Reynolds  number.  For 
an  optimized  cavity  length  for  the  pressure  variation  case,  this 

g 

Reynolds  number  is  in  the  neighborhood  of  30  x 10  . 

It  was  found  that  the  rate  of  pressure  recovery  is  different 
from  the  rate  of  velocity  profile  recovery.  The  pressure  recovers 
faster  (that  is  with  less  relaxation  distance)  than  does  the 
velocity  profile.  The  reason  for  these  differences  is  that  the 
surface  pressure  is  mainly  determined  by  the  immediate  external 
flow  field  while  the  change  in  the  velocity  profile  depends  on 
the  mass  entrainment  from  the  external  flow.  The  mass  entrain- 
ment through  mixing  during  the  rehabilitation  process  requires  a 
longer  relaxation  distance.  At  an  optimized  cavity  length  for 
the  pressure  recovery,  it  was  found  that  the  final  peak  pressure 
as  well  as  the  reseparation  pressure  were  higher  than  those  for 
the  no  upstream  disturbance  flow  case.  At  an  optimized  cavity 
length  for  the  velocity  profile  recovery,  it  was  found  that  the 
skin  friction  can  be  higher  than  that  for  the  case  with  no  up- 
stream disturbance. 

By  observing  the  measured  velocity  profile  development 
during  the  rehabilitation  phase,  it  was  found  that  the  recovery 
(the  pressure  as  well  as  the  velocity)  starts  from  the  wall  and 
spreads  outward.  The  inner  portion  of  the  outer  boundary  layer 
(i.e.,  the  wake- component)  has  a longer  memory  of  the  upstream 
disturbance  than  does  the  outer  layer  (i.e.,  the  wall-component). 
These  points  are  even  more  clearly  demonstrated  by  analyzing 
the  velocity  defect  and  the  law  of  wall  and  wake. 

It  was  found  that  the  value  used  for  the  law  of  the  wall 
in  the  flow  over  a flat  plate  does  not  agree  with  the  condition 
where  a large  pressure  gradient  is  presented.  Therefore,  the 
slope  on  Coles'  plot  deviates  from  the  universally  accepted 
value  near  the  reattachment  and  the  reseparation  points . 
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A method  of  computing  the  shape  of  the  inviscid  streamlines 
has  been  developed  for  subsonic  flows  past  cavities  using  the 
experimentally  measured  data.  Using  such  inviscid  streamlines, 
inviscid  transonic  flow  computation  for  flows  past  cavities  at 
transonic  speeds  yields  pressure  distribution  on  the  cavities 
which  agrees  with  the  experimentally  measured  pressure  distri- 
butions at  transonic  speeds.  This  indicates  that  the  extent  of 
the  viscous  region  does  not  change  much  as  the  Mach  number  is 
increased  from  subsonic  to  transonic  level. 

A new  finite  difference  technique  for  solving  the  viscous- 
transonic  equations  has  also  been  developed  which  includes  the 
effects  of  compressive  viscosity  in  the  regions  of  shocks.  The 
method  has  been  applied  for  solving  flows  past  symmetric  airfoils 
and  it  has  been  established  that  the  results  agree  with  those 
obtained  by  other  methods  and  the  present  technique  has  certain 
computational  advantages.  Using  this  method,  pressure  distri- 
butions on  circular  airfoils  with  a pre-assigned  separation  wake 
either  just  behind  the  point  where  the  shock  touches  the  airfoil 
or  at  a point  further  downstream  have  been  computed. 

This  work  has  gone  some  way  towards  developing  an  under- 
standing of  transonic  strong  interaction  flows  but  significant 
additional  work  is  still  required  for  a complete  knowledge  of  the 
problem. 
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APPENDIX 

NUMERICAL  SOLUTION  OP  VISCOUS-TRANSONIC  EQUATION 

The  small  disturbance  equations  for  the  disturbance 
velocity  potential  for  inviscid  transonic  flow  past  an  airfoil 
with  effects  of  compressive  viscosity  included  are 
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where  is  the  inviscid  transonic  similarity  parameter  and 
the  viscous  transonic  similarity  parameter  defined  by 
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6 is  the  thickness  to  chord  ratio  of  the  airfoil  and  Re  and  Pr 
are  the  Reynolds  number  and  Prandtl  number  respectively,  based 
on  longitudinal  viscosity  and  half-chord  length.  x,y  and  are 
suitably  non-dimens ionalized  variables. 

Boundary  conditions  are 
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Introducing  U - *x  and  W = *xx>  Equation  (1)  is  written  as  a 
system  of  three  equations  which  have  only  first  order  derivatives 

T T*1  V 


m x. 


+x  - u 

Ux  = W 


(Kvw  + Ktn  - %u2)x  - -+ 


The  third  equation  is  written  in  the  divergence  form  so  that  the 

corresponding  difference  equations  would  preserve  the  conservation 
laws . 

Difference  Notation: 


Xi  “ xi-l  + ht  ’ i = 2,3,  , . , I 


x i “ -«  , xt  - +03 
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Yj  = y^-i+  Sj  . j = 2,3,  . . , J 

= 0,  yx  = gj/2  , ya  = + » 

ui,j  - 


i-1  i 


The  difference  equations  corresponding  to  the  equations 
(6)  are  written  at  the  point  The  difference  scheme 

is  second  order  accurate  and  the  grid  spacing  is  non-uniform  in 
both  x and  y directions.  The  difference  scheme  is  given  by 
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equation  (4) . 

This  difference  scheme  which  is  similar  to  the  "box  scheme" 
developed  by  Keller  (Ref.  32)  for  parabolic  problems,  would 
remain  the  same  in  both  subsonic  and  supersonic  regions  of  the 
flow  field,  unlike  the  mixed  difference  schemes  used  by  most  others 
for  inviscid  transonic  flow  computations.  The  equations  (7) 
together  with  the  boundary  conditions  are  a set  of  non-linear 
algebraic  equations  for  <j> , U and  W at  all  the  mesh  points  (i,j); 

I;j  = l,2,  . ..  J.  These  equations  are  solved  by  a 
horizontal  line  relaxation  method.  Assuming  some  approximation 
for  the  solution  in  the  entire  region,  equations  on  line  j = 1 
are  solved  for  obtaining  the  next  approximation  for  the  variables 
on  that  line.  However,  even  on  this  line,  equations  are  non- 
linear and  they  themselves  are  solved  by  Newton's  iteration. 

T7hile  carrying  out  those  Newton  iteration  steps,  structure  of 
the  equation  matrix  which  is  a tri-diagonal  block  matrix  with 
additional  zeros  in  some  of  the  blocks,  is  taken  advantage  of. 

Then  equations  on  line  j = 2 are  solved  and  so  on.  One  sweep  is 
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completed  after  the  equations  on  all  the  horizontal  lines  are 
solved.  This  gives  next  approximation  in  the  entire  field  and 
the  next  sweep  is  carried  out.  It  takes  several  hundred  sweeps 
for  the  relaxation  scheme  to  converge  if  the  starting  approxi- 
mation is  very  far  from  the  actual  solution.  However,  when  one 
solution  is  obtained  for  certain  values  of  the  parameters  of  the 
problem,  that  solution  is  used  as  the  starting  approximation  for 
obtaining  solutions  for  slightly  different  values  of  the  para- 
meters. Convergence  would  be  rapid  in  such  cases. 

In  a typical  run,  about  125  points  in  the  x-direction  and 
about  30  points  in  the  y-direction  are  used.  The  mesh  spacing 
is  non-uniform  and  mesh  lines  are  introduced  in  the  field  when 
necessary  in  regions  of  strong  shocks  for  better  resolution  during 
the  computation,  and  also  lines  can  be  taken  out  of  computation 
when  it  is  desirable. 
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C£ 

c 

p 

f (x) 


r (x) 


Nomenclature 

Sound  speed 

Wall-law  intercept 

Skin  friction  coefficient 

Pressure  coefficient 

Cavity  or  airfoil  shape  function 

Clauser's  velocity  defect  parameter 

Step  height 

Transonic  similarity  parameter 

Inviscid  transonic  similarity  parameter 

Viscous  transonic  similarity  parameter 

Cavity  length 

Reattachment  distance 

Separation  distance 

Distance  to  forward  facing  step 

Distance  to  rearward  facing  step 

Mach  number 

Pressure 

Prandtl  number 

Recovery  factor 

{(l-x)/(l+x)}% 

Reynolds  number 

Axial  velocity 

Friction  velocity 

van  Driest 's  generalized  velocity 

Axial  distance 

Distance  downstream  from  rearward  facing  step 
Distance  upstream  from  forward  facing  step 
Transverse  coordinate 

1/3 

Stretched  y coordinate,  6 y 
Angle  of  attack 
Specific  heat  ratio 

Boundary  layer  thickness;  cavity  depth  or  airfoil 
thickness  to  chord  ratio. 
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6* 

6* 
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n 

A 

e 

e 
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K 

P 

V 

c 

TT 
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to 


Boundary  layer  displacement  thickness 
Cole's  displacement  thickness 
Integration  variable 
Clauser's  thickness  parameter 
Boundary  layer  momentum  thickness 
Cole’s  momentum  thickness 
von  Karman  constant 
Coefficient  of  viscosity 
Coefficient  of  kinematic  viscosity 
Integration  variable 
Cole's  profile  parameter 
Density 

Disturbance  velocity  potential 
Wake  function 


SUBSCRIPTS : 
b Base 

e Edge  of  boundary  layer 

i Characteristic 

Distance  to  forward  facing  step 
Lr  Distance  to  rearward  facing  step 

i/2  Center  of  cavity 

m Peak 

p Plateau 

r Reattachment 

ref  Reference 

s Separation 

t Total 

w Wall 

x Axial  distance 

6 Boundary  layer  thickness 

00  Free  stream 
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